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SUMMARY 


The  results  of  a  continuing  program  of  theory  and  experiment 
concerned  with  the  aerodynamic  characteristics  of  airfoils  in 
nonuniformly  sheared  flows  are  presented.  Development  of  a  digital 
computer  program  to  calculate  airfoil  pressure  distributions  in  a 
particular  type  of  two-dimensional,  inviscid  flow  with  nonuniform 
shear  was  completed,  and  airfoil  pressure  distributions  were  computed 
for  specific  configurations. 

Wind-tunnel  tests  were  performed  with  an  airfoil  in  a  two- 
dimensional,  nonuniform  flow  with  approximately  the  same  free-stream 
velocity  profile  as  that  used  in  the  computer  program.  Airfoil  pressure 
distributions  were  obtained  with  the  airfoil  at  several  locations  in  the  non- 
uniform  flow,  and  throughout  a  range  of  angles  of  attack  including  angles 
at  which  complete  separation  occurred.  It  was  found  that  the  separation 
characteristics  varied  with  the  location  of  the  airfoil  in  the  sheared 
flow.  Also,  promotion  or  delay  of  separation  was  related  qualitatively 
to  variations  in  the  adverse  pressure  gradient  on  the  airfoil's  upper 
surface;  and  it  was  observed  that  the  trends  are  similar  to  those  found 
in  uniform  flow.  Finally,  it  was  found  that  if  lift  coefficients  referenced 
to  stagnation  streamline  free-stream  dynamic  pressures  \_(^i  )s  ^are 
plotted  against  angles  of  attack  measured  from  zero  lift  (at),  the  resulting 
data  fall  into  a  relatively  narrow  band  centered  on  the  uniform  flow  experi¬ 
mental  cc  versus  at  curve. 

Computed  pressure  distributions  were  in  good  agreement  with 
the  corresponding  experimentally  measured  pressure  distributions.  It 
would  appear,  therefore,  that  observed  anomalous  aerodynamic  behavior 
near  flow  separation  of  airfoils  in  nonuniform  flows  is  a  direct  result  of 
the  effects  of  the  inviscid  sheared  flow  on  the  airfoil  pressure  distribution. 


The  validated  analytical  technique  for  computing  pressure  distributions 
on  an  airfoil  in  an  inviscid  nonuniform  flow  provides  a  powerful  tool 
for  investigation  of  such  anomalous  behavior. 
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Cornell  Aeronautical  Laboratory,  Inc.  (CAL),  Buffalo,  New  York,  for 
the  U.  S.  Army  Aviation  Materiel  Laboratories  (USAAVLABS),  Fort 
Eustis,  Virginia.  The  work  was  performed  during  the  12-month  period 
starting  in  January  1965  under  Contract  DA  44-177-AMC-268(T), 
"Theory  of  Airfoils  in  Nonuniform  Sheared  Flows".  Mr.  P.  Cancro 
administered  the  project  for  USAAVLABS.  This  work  is  a  continuation 
of  a  research  program  carried  out  at  CAL  over  a  period  of  years  from 
1959  under  Contracts  DA  44-177-TC-439  and  DA  44-177-AMC-70(T). 
Mr.  W.  G.  Brady  and  Dr.  G.  R.  Ludwig  are  authors  of  this  report. 

Mr.  Brady's  primary  responsibility  was  the  analytical  work,  and 
Dr.  Ludwig's  primary  responsibility  was  the  experimental  program. 
The  invaluable  contributions  of  Mr.  J.  Nemeth  and  Mr.  C.  Ryan  of 
CAL  are  appreciated. 


v 


BLANK  PAGE 


CONTENTS 


Page 

SUMMARY  iii 

FOREWORD  v 

LIST  OF  ILLUSTRATIONS  viii 

LIST  OF  SYMBOLS  xii 

INTRODUCTION  1 

THEORETICAL  PROGRAM  4 

EXPERIMENTAL  PROGRAM  14 

DISCUSSION  OF  EXPERIMENTAL  RESULTS  21 

COMPUTED  RESULTS  AND  COMPARISON  WITH 

EXPERIMENTAL  RESULTS  36 

CONCLUSIONS  AND  RECOMMENDATIONS  44 

REFERENCES  88 

DISTRIBUTION  90 

APPENDIX.  MATHEMATICAL  FORMULATION  AND 
IMPLEMENTATION  OF  SHEARED  FLOW 

AERODYNAMICS  COMPUTER  PROGRAM  91 


ILLUSTRATIONS 


Figure  Page 

1  Examples  of  Free-Stream  Velocity  Profiles 

Represented  by  Piecewise  Linear  Segments.  47 

2  Schematic  of  Flow  Model  for  IBM  Computer  Program.  47 

3  Pressure  Coefficient,  Cp  ,  Versus  Fraction  of  Chord, 

%/c  ,  for  Uniform  Flow  From  Potential  Flow  Theory, 
and  From  Computer  Program  Without  Smoothing; 

os  =  10.  0  Degrees.  48 

4  Pressure  Coefficient,  Cp  ,  Versus  Fraction  of  Chord, 

X/c ,  for  Uniform  Flow  From  Potential  Flow  Theory, 
and  From  Computer  Program  With  Smoothing 

(  R  r  0.  10  );  cc  -  10.  0  Degrees.  49 

5  View  of  Airfoil  Model  Installed  in  Wind-Tunnel  Test 

Section.  50 

6  Shear  Screen  Assembly.  51 

7  Measured  Free-Stream  Velocity  Distribution  in  Non- 
uniformly  Sheared  Flow  at  Model  Midchord  Station.  52 

8  Measured  Free-Stream  Flow  Angularity  in  Nonuni- 

formly  Sheared  Flow  at  Model  Midchord  Station.  53 

9  Section  Lift  Coefficient,  ,  Versus  Angle  of  Attack, 

OC,  for  Uniform  Flow.  54 

10  Pressure  Coefficient,  Cp  ,  Versus  Fraction  of  Chord, 

x>/c  ,  Uniform  Flow,  (X  =  0  Degrees.  55 

11  Pressure  Coefficient,  Cp ,  Versus  Fraction  of  Chord, 

x,/c  *  Uniform  Flow,  (X  -  7.  5  Degrees.  56 

12  Pressure  Coefficient,  CF ,  Versus  Fraction  of  Chord, 

%/c  *  Uniform  Flow,  CC  =  10.  0  Degrees.  57 


vm 


Figure 


Page 


13  Pressure  Coefficient,  Cp,  Versus  Fraction  of  Chord, 

X/c ,  Uniform  Flow,  Low  Solidity  Turbulence  Screen, 

06  =  13.0  Degrees,  U0  =  66  Ft.  /Sec.  58 

14  Pressure  Coefficient,  Cp  ,  Versus  Fraction  of  Chord, 
tC/c,  Uniform  Flow,  Low  Solidity  Turbulence  Screen, 

<X  -  16.  5  Degrees,  U0  -  66  Ft.  /Sec.  59 

15  Pressure  Coefficient,  Cp  ,  Versus  Fraction  of  Chord, 

%/c ,  Uniform  Flow,  Low  Solidity  Turbulence  Screen, 

OS  *  16.  7  Degrees,  0o  =  66  Ft.  /Sec.  60 

16  Section  Quarterchord  Moment  Coefficient,  cmcH  ,  Versus 

Angle  of  Attack,  a, in  Uniform  Flow  With  Low  Solidity 
Turbulence  Screen,  61 

17  Section  Pressure  Drag  Coefficient,^, Versus  Angle  of 

Attack,  cct  in  Uniform  Flow  With  Low  Solidity  Turbulence 
Screen.  62 

18  Nondimensional  Section  Lift,  Cj,  Versus  Angle  of  Attack, Or, 

Nonuniformly  Sheared  Flow.  63 

19  Section  Lift  Coefficient  Referenced  to  Midchord  Free- 
Stream  Dynamic  Pressure,  (Ojip,  Versus  Angle  of 

Attack,  cc  ,  Nonuniformly  Sheared  Flow.  65 

20  Section  Lift  Coefficient  Referenced  to  Midchord  Free- 
Stream  Dynamic  Pressure,  ( cj)M ,  Versus  Angle  of 

Attack,  06  ,  Nonuniformly  Sheared  Flow  (Reference  6).  67 

21  Section  Lift  Coefficient,  (c^)^.  Versus  Section  Pressure 

Drag  Coefficient,  (e^)*,  Nonuniformly  Sheared  Flow.  69 

22  Section  Quarterchord  Moment  Coefficient,  (cmeM  )*  , 

Versus  Section  Lift  Coefficient,  (cj  )*,  Nonuniformly 
Sheared  Flow.  71 


23  Pressure  Coefficient  Referenced  to  Midchord  Free- 

Stream  Dynamic  Pressure,  (C?)Kt  Versus  Fraction  of 
Chord,  x/c  ,  for  Uniform  Flow  and  Nonuniformly 
Sheared  Flow. 


73 


Figure 


24(a)  Section  Lift  Coefficient  Referenced  to  Free-Stream 
Dynamic  Pressure  of  Stagnation  Streamline, (Cj)s, 
Versus  Angle  of  Attack  Measured  From  Zero  Lift, 
Nonuniformly  Sheared  Flow, 


Page 


74 


24(b)  Section  Lift  Coefficient  Referenced  to  Midchord  Free- 
Stream  Dynamic  Pressure,  [cj  )*,  Versus  Angle  of 
Attack  Measured  From  Zero  Lift,  oc ,  Nonuniformly 
Sheared  Flow,  75 

25  Pressure  Coefficient  Referenced  to  Free-Stream 
Dynamic  Pressure  of  Stagnation  Streamline,  (Cp)5% 

Versus  Fraction  of  Chord, x/c,  for  dc  -  5  Degrees, 

Uniform  Flow  and  Nonuniformly  Sheared  Flow,  76 

26  Pressure  Coefficient  Referenced  to  Free-Stream 
Dynamic  Pressure  of  Stagnation  Streamline,  (Cp)% , 

Versus  Fraction  of  Chord ,x/c,  for  a  =  12.5  Degrees, 

Uniform  Flow  and  Nonuniformly  Sheared  Flow.  77 

27  Displacement  of  Nonuniformly  Sheared  Flow  Airfoil 
Stagnation  Streamline  Relative  to  Airfoil  Midchord, 

h  =  -  3  inches.  78 

28  Displacement  of  Nonuniformly  Sheared  Flow  Airfoil 
Stagnation  Streamline  Relative  to  Airfoil  Midchord, 

h  =  -  2  inches.  79 

29  Displacement  of  Nonuniformly  Sheared  Flow  Airfoil 
Stagnation  Streamline  Relative  to  Airfoil  Midchord, 

h  -  -  1  inch.  80 

30  Displacement  of  Nonuniformly  Sheared  Flow  Airfoil 
Stagnation  Streamline  Relative  to  Airfoil  Midchord, 

h  =  +  1  inch.  81 

31  Displacement  of  Nonuniformly  Sheared  Flow  Airfoil 
Stagnation  Streamline  Relative  to  Airfoil  Midchord, 

h  =  +  3  inches.  82 

32  Comparison  of  Computed  and  Measured  Pressure 
Coefficients,^*)*,  Versus  Fraction  of  Chord, Z/c  , 

Nonuniformly  Sheared  Flow,  h  =  -  2  inches,  <x>  - 

7.  1  Degrees.  83 


x 


Figure 


Page 

33  Comparison  of  Computed  and  Measured  Pressure 
Coefficients,^^  ,  Versus  Fraction  of  Chord,  x/c , 
Nonuniformly  Sheared  Flow,  h  =  -  2  inches,  Cc 

9.  5  Degrees.  84 

34  Comparison  of  Computed  and  Measured  Pressure 
Coefficients,  (Cp)g,  Versus  Fraction  of  Chord,  x/c  , 
Nonuniformly  Sheared  Flow,  h  =  +  3  inches,  oc  « 

8.  2  Degrees.  35 

35  Comparison  of  Computed  and  Measured  Section 
Lift  Coefficient  Referenced  to  Midchord  Free- 
Stream  Dynamic  Pressure,  Nonuniformly  Sheared 

Flow.  86 

36  Comparison  of  Computed  and  Measured  Section 
Quarterchord  Moment  Coefficient,  Nonuniformly 

Sheared  Flow.  87 


xi 


LIST  OF  SYMBOLS 


c 


<>d 

(c*)* 

ci 

<cj)* 

(tj)s 


Zj 


W4 


(Cme/v)/f 

CP 

(V* 


(C,)s 


4 

4r 

m 

i(*> 


Airfoil  chord,  inches. 

Uniform  flow  section  drag  coefficient;  Cj  ®  — : — — 77 - - 

6  jjO(0J*c 

Sheared  flow  section  drag  coefficient,  referred  to  free-stream 
velocity  at  airfoil  midchord  position;  (c^)*  =  UN,r  sFA*  • 

Uniform  flow  section  lift  coefficient;  c»  -  UF?  • 

Sheared  flow  section  lift  coefficient,  referred  to  free-stream 
velocity  at  airfoil  midchord  position;  *c' " ’ 

Sheared  flow  section  lift  coefficient,  referred  to  free-stream^^ 
velocity  along  airfoil  stagnation  streamline;  (cj)a  *  -  - - 

Sheared  flow  section  lift  coefficient,  referred  to  average  free- 

stream  velocity  in  wind-tunnel  section;  Cj  =  ..yyjL*.?**  ■ 

y  J  iP(U)‘c 

Uniform  flow  section  moment  coefficient  about  the  airfoil 
quarter  chord,  positive  nose  up;  c„ 


moment  pe/t  uNtr  spam 


'c/4 


Sheared  flow  section  moment  coefficient  about  the  airfoil 
quarterchord,  positive  nose  up;  (c„c/^y  *±**y^ic* 

Uniform  flow  pressure  coefficient;  C„  «•  —  . 

p 

Sheared  flow  pressure  coefficient,  referred  to  free-stream 
velocity  at  airfoil  midchord  position;  (C^R  -  ‘ 

Sheared  flow  pressure  coefficient,  referred  to  free-stream 

--  r~Po  . 


velocity  along  airfoil  stagnation  streamline;  (Cp)s  f  (U  y- 

Distance  upstream  and  downstream  of  airfoil  midchord  where 
sheared  flow  dividing  streamline  is  assumed  to  be  at  its  free- 
stream  location  (see  Figure  2),  inches. 


Distance  upstream  and  downstream  of  airfoil  midchord  where 
wind-tunnel  wall  vortex  sheet  singularity  strengths  are  assumed 
to  have  reached  their  free-stream  values  (see  Figure  2),  inches. 

Function  defined  by  Equation  (1). 

Displacement  of  dividing  streamline  (see  Figure  2),  inches. 


Xll 


A  Vertical  height  of  airfoil  midchord  above  tunnel  centerline,  inches. 
K  Shear  parameter;  K  *  * 

/T  Smoothing  parameter;  see  Equation  (A-44). 

L  Vertical  distance  from  tunnel  centerline  to  wall  (see  Figure  2), 
inches. 


p(*.y)  Static  pressure  at  (x,y  ),  pounds  per  square  foot. 

p0  Static  pressure  in  free  stream,  pounds  per  square  foot. 

U 9  Uniform  flow  free-stream  velocity,  feet  per  second. 

U(y)  Sheared  flow  free-stream  velocity,  feet  per  second. 

U%  Sheared  flow  free-stream  velocity  at  y  =  h  ,  feet  per  second. 

Us  Sheared  flow  free-stream  velocity  of  stagnation  streamline, 

feet  per  second. 

Uf  Free-stream  velocity  at  lower  wind-tunnel  wall,  feet  per  second. 

Uz  Free-stream  velocity  at  upper  wind-tunnel  wall,  feet  per  second. 

Free-stream  velocity  along  vertical  centerline  in  wind  tunnel, 
feet  per  second. 

U  Average  free-stream  velocity  in  wind-tunnel  test  section,  feet 
per  second. 

u>  Velocity  component  in  X, -direction,  feet  per  second. 

V  Velocity  component  in  y -direction,  feet  per  second. 

X  Coordinate  along  axis  parallel  to  wind-tunnel  walls,  positive 
downstream;  x>  =  0  at  airfoil  midchord,  inches. 

tL  ^-coordinate  of  airfoil  leading  edge,  inches. 

^-coordinate  of  airfoil  trailing  edge,  inches. 

AYS  Vertical  distance  above  the  airfoil  midchord  reference  height, 
h,  where  U(y)  -  Us  in  the  free-stream  velocity  profile,  inches, 
(see  sketch,  p.  33). 

y  Coordinate  perpendicular  to  x-axis,  positive  up;  <J  -  0  at  wind- 
tunnel  centerline  ,  inches. 


xiii 


Airfoil  angle  of  attack  measured  from  zero  lift  angle  of  attack, 
positive  leading  edge  up,  degrees. 

Airfoil  angle  of  attack,  positive  leading  edge  up,  degrees. 

Singularity  strength  of  vortex  sheet  representing  airfoil 
aft  lower  surface. 

Singularity  strength  of  vortex  sheet  representing  airfoil 
aft  upper  surface. 

Singularity  strength  of  vortex  sheet  representing  airfoil 
leading  edge  surface. 

Singularity  strength  of  vortex  sheet  representing  lower 
wind-tunnel  wall. 

Singularity  strength  of  vortex  sheet  representing  upper 
wind-tunnel  wall. 

Air  density,  assumed  constant,  slugs  per  cubic  foot. 

Stream  function;  see  Equation  (1). 


xiv 


INTRODUCTION 


Many  of  the  analytical  problems  of  wing  aerodynamics  in  flows 
at  low  subsonic  speeds  are  treated  successfully  within  the  framework  of 
potential  flow  theory.  However,  potential  flow  theory  is  inadequate  for 
predicting  analytically  the  interaction  between  a  wing  and  a  propeller 
slipstream  when  the  slipstream  intersects  the  wing  and  the  flow  in  the 
propeller  slipstream  is  not  irrotational,  i.  e.  ,  when  the  slipstream  flow 
has  significant  shear.  An  additional  difficulty  which  gives  rise  to 
nonlinearities  in  the  analytical  treatment  results  from  the  interaction  of 
the  free  boundaries  of  the  slipstream  flow  and  the  wing. 

A  number  of  analytical  treatments  of  the  aerodynamics  of  sheared 
flows  have  been  published;  typical  are  References  1  through  4.  Exact 
aerodynamic  theories  for  uniform  two-dimensional  shear  (constant 
vorticity)  are  available  (References  1  and  4),  but  the  available  solutions 
for  nonuniform  shear  (Reference  3,  for  example)  are  not  particularly 
useful  to  the  aerodynamicist  because  of  limitations  of  the  theoretical 
assumptions.  These  are  primarily  restrictions  to  small  shear  and 
shear  gradient.  Three-dimensional  sheared  flow  theory  (Reference  2, 
for  example)  has  hardly  advanced  beyond  attempts  to  arrive  at  a  formu¬ 
lation  which  is  mathematically  tractable. 

In  view  of  this,  the  aerodynamicist  is  confronted  with  a  difficult 
task  if  rational  design  data  are  to  be  provided  for  the  low-speed  aero¬ 
dynamic  characteristics  of  VTOL  or  STOL  configurations  in  which  much 
of  the  wing  is  immersed  in  a  propeller  wake.  Reliance  must  be  placed 
on  experimental  data,  either  from  wind-tunnel  tests  or  from  free -flight 
data  from  prototype  aircraft.  In  either  case,  without  the  foundation  of 
prior  experience  and  theoretical  understanding,  much  guesswork  and  many 
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painful  and  expensive  attempts  at  c.ut-and-try  fixes  to  problems  in 
performance,  stability,  and  handling  characteristics  may  be  required 
to  produce  a  satisfactory  configuration. 

The  research  reported  here  is  part  of  a  continuing  program 
concerned  with  the  aerodynamic  characteristics  of  sheared  flows. 

Initial  results  were  obtained  during  a  previous  program  of  theoretical 
and  experimental  research  on  low-speed  aerodynamics  relevant  to 
STOL  and  VTOL  aircraft  conducted  for  the  U.  S.  Army  Aviation  Materiel 
Laboratories  (formerly  USATRECOM).  References  4,  5,  6,  and  7  report 
the  results  of  that  part  of  the  study  devoted  to  sheared  flows.  Of 
particular  interest  was  the  discovery  that  nonuniform  shear  in  a  two- 
dimensional  flow  could  have  a  marked  effect  on  airfoil  stall  character¬ 
istics  (Reference  6).  It  is  noted  in  Reference  6  that  in  the  specific  two- 
dimensional  nonuniformly  sheared  flow  treated,  large  variations  in 
maximum  lift  were  obtained  near  the  flow  centerline  when  the  airfoil 
vertical  position  was  varied  by  distances  of  the  order  of  magnitude  of 
the  airfoil  thickness. 

The  work  reported  in  Reference  8  continued  the  preliminary 
experimental  program  of  Reference  7, which  was  concerned  with  the 
effects  of  shear  on  a  two-dimensional  airfoil  immersed  in  a  three- 
dimensional,  axially-symmetr ic ,  nonuniformly  sheared  flow.  Also 
reported  in  Reference  8  was  the  initial  effort  on  development  of  an 
analytical  technique  for  computing  the  aerodynamic  characteristics  of 
an  airfoil  in  two-dimensional,  inviscid,  nonuniformly  sheared  flow. 

The  research  reported  here  continued  the  development  of  this 
analytical  technique.  A  digital  computer  program  was  developed  to 
compute  pressure  distributions  on  a  two-dimensional  airfoil  in  a 
particular  nonuniformly  sheared  flow.  Experimental  airfoil  pressure 
distribution  data  were  obtained  in  the  wind  tunnel  for  a  two-dimensional, 


2 


free-stream  velocity  profile  nearly  the  same  as  that  in  the  theory. 

These  data  were  obtained  to  enable  the  accuracy  of  the  analytical 
techniques  to  be  verified.  It  was  expected  that  the  availability  of  an 
accurate  theory  would  be  of  great  value  in  investigations  to  determine 
the  mechanism  responsible  for  the  unusual  stall  characteristics  of  the 
airfoil  in  nonuniformly  sheared  flows,  as  well  as  to  obtain  a  more  general 
understanding  of  sheared  flow  aerodynamic  phenomena. 
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THEORETICAL  PROGRAM 


BACKGROUND  -  FUNDAMENTAL  THEORETICAL  CONSIDERATIONS 

The  research  program,  of  which  the  work  reported  here  is  a 
part,  is  motivated  by  the  lack  of  a  general  understanding  of  sheared 
flow  aerodynamics.  In  particular,  an  understanding  of  the  mechanisms 
inherent  in  the  maximum  lift  behavior  of  airfoils  in  sheared  flows 
should  enable  the  aircraft  designer  to  take  advantage  of  this  behavior. 

For  example,  experimental  results  obtained  in  Reference  6  indicate 
that  large  variations  in  maximum  lift  can  occur  with  small  changes  in 
vertical  position  on  an  airfoil  in  two-dimensional,  nonuniformly  sheared 
flows  with  large  shear  and  shear  gradients. 

The  unusual  and  interesting  maximum  lift  behavior  in  a  particular 
nonuniformly  sheared  flow,  demonstrated  experimentally  in  Reference  6, 
could  arise  from  any  one  or  a  combination  of  a  number  of  possible 
mechanisms.  Among  those  possibilities  are  the  following: 

1.  Subtle  wind-tunnel  wall  effects  associated  with  the  sheared 
flow. 

2.  An  altering  of  the  airfoil  pressure  distribution,  solHv  on 
the  basis  of  the  inviscid  sheared  flow,  so  as  to  modify  the 
separation  of  the  airfoil  boundary  layer. 

3.  An  interaction  of  the  boundary  layer  with  the  inviscid 
sheared  flow. 

It  might  seem  reasonable  to  test  the  various  hypotheses  selectively 
for  their  relative  importance  by  means  of  carefully  controlled  wind- 
tunnel  experiments;  however,  this  approach  can  be  time-consuming 
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and  expensive.  The  large  number  of  potentially  important  variables 
and  the  probable  importance  of  nonlinearities  make  an  all-experimental 
type  of  investigation  even  less  attractive.  If  an  accurate  inviscid  aero¬ 
dynamic  theory  for  nonuniformly  sheared  flow  were  available,  however, 
the  experimenter  would  have  a  very  powerful  added  tool  to  use  in  con¬ 
junction  with  carefully  chosen  experiments.  In  particular,  such  an 
inviscid  theory  should  enable  the  experimenter  to  identify,  or  at  least 
to  isolate,  viscous  effects.  If,  at  the  same  time,  the  theory  is  capable 
of  predicting  wind-tunnel  wall  effects,  then  an  additional  uncertainty  in 
the  test  program  could  be  eliminated. 

The  theoretical  study  initiated  during  the  work  reported  in 
Reference  8,  and  continued  during  the  present  program,  has  the 
ultimate  objective  of  providing  an  analytical  technique  for  predicting 
aerodynamic  characteristics  of  airfoils  in  arbitrary  two-dimensional, 
inviscid,  nonuniformly  sheared  flows.  The  assumptions  and  conditions 
underlying  the  investigation  are  as  follows: 

1.  The  flow  and  the  airfoil  section  are  two-dimensional. 

2.  The  flow  is  inviscid  and  incompressible. 

3.  The  mathematical  flow  model  must  allow  for  large  shear  and 
shear  gradients  in  the  free  stream. 

4.  The  mathematical  flow  model  must  be  consistent  with  large 
disturbances  in  the  flow  field. 

The  reason  for  specifying  condition  4  above  is  related  to  the  experi¬ 
mental  findings  that  a  sheared  flow  can  either  promote  or  retard  flow 
separation  at  large  angles  of  attack.  Therefore,  the  large  disturbances 
corresponding  to  these  large  angles  of  attack  must  be  tolerated  by  the 
mathematical  model.  The  linearizing  techniques  used  in  thin  airfoil 
theory,  therefore,  would  not  be  valid  for  those  cases  in  which  high 
angles  of  attack  are  considered. 


The  main  characteristic  of  a  sheared  flow  which  distinguishes 
it  from  a  nonsheared  flow  is  the  fact  that  it  is  rotational.  The  powerful 
techniques  of  potential  theory,  including  conformal  mapping,  which  can 
be  used  for  two-dimensional, ir rotational  flow  problems,  are  generally 
not  applicable  to  rotational  flows. 

Rather  than  Laplace's  equation,  the  stream  function  ^  must 

satisfy 

-  f(*)  \  (1) 

i.  e.  ,  V*y  is  a  constant  along  a  streamline  (along  which  the  stream 
function  is  also  a  constant).  It  can  be  shown  that,  in  two-dimensional 
flow,  this  constant  is  the  vorticity;  i.  e. ,  the  vorticity  is  a  constant 
along  a  streamline  in  an  inviscid  two-dimensional,  rotational  flow. 
Except  for  the  case  of  uniform  shear  ( £( f)  -  constant),  or  the 
exceptional  case  f(f)  =  cp ,  c  =  constant,  one  is  faced  with  solving 
a  nonlinear,  partial  differential  equation. 

A  general  survey  of  available  theory  for  the  aerodynamics  of 
sheared  flows  was  presented  in  Reference  8.  It  was  concluded  there 
that  the  available  theories  were  inadequate  for  the  purposes  of  the 
present  research.  In  particular,  the  requirements  that  a  suitable 
theory  permit  large  shear  and  shear  gradients  and  be  consistent  with 
large  disturbances  in  the  flow  field  were  not  satisfied.  In  view  of  this, 
a  new  approach  was  adopted,  based  on  a  numerical  technique.  Although 
the  resultant  digital  computer  program  was  outlined  in  Reference  8, 
this  outline  is  repeated  here  for  the  sake  of  completeness.  The  basic 
ideas  behind  the  implementation  of  the  computer  program  are  presented 
in  the  body  of  this  report.  The  fundamental  equations  and  details  of 
the  actual  program  which  was  developed  for  use  on  an  IBM  7044  digital 
computer  are  presented  in  the  appendix  of  this  report. 
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A  digital  computer  program  for  computing  an  axially  symmetric, 
free-streamline  flow  (finite  impinging  jet)  was  reported  in  Reference  9. 

The  success  of  this  program  indicated  a  method  of  approach  which  should 
permit  the  calculation  of  the  aerodynamic  characteristics,  including 
pressure  distributions,  of  any  two-dimensional  airfoil  at  arbitrary  angle 
of  attack  (in  nonseparated  flow)  in  two-dimensional,  nonuniformly  sheared 
flows  with  arbitrary  shear.  The  only  restriction  appears  to  be  that  the 
free-stream  velocity  profile  is  such  that  it  can  be  approximated  by  piece - 
wise  linear  segments.  Examples  of  such  profiles  are  shown  in  Figure  1. 

It  should  be  recalled  that  the  fundamental  difficulty  in  treating 
two-dimensional,  inviscid,  nonuniformly  sheared  flows  analytically  is 
related  to  the  fact  that  the  fundamental  governing  equation  (Equation  (1)) 
is  nonlinear  for  general  nonuniformly  sheared  flows.  This  equation 
states  that  the  vorticity  in  the  flow  is  a  constant  along  streamlines; 
inasmuch  as  the  shape  of  the  streamline  is  one  of  the  unknowns  in  the 
problem,  the  functional  form  of  the  right-hand  side  of  Equation  (1)  is  not 
known  beforehand.  If  the  free-stream  velocity  profile  can  be  approximated 
by  piecewise  linear  segments,  as  in  Figure  1,  then  the  streamlines 
passing  through  those  points  at  which  the  velocity  gradient  changes 
(points  (a),  (b),  and  (c)  in  Figure  1,  for  example)  separate  regions  of 
constant  vorticity  throughout  the  flow.  Thus,  the  relationship  to  the  free- 
streamline  flow  problem  (free -boundary  problem)  is  clear;  as  part  of  the 
solution  of  the  problem,  the  shape  of  thc«»e  streamlines  separating 
regions  of  constant  vorticity  must  be  determined.  In  effect,  a  problem 
with  a  nonlinear  governing  equation  and  known  boundary  conditions  is 
replaced  by  a  problem  with  a  linear  governing  equation  (  V*\f/  =  constant) 
within  regions  with  free  (and,  hence,  unknown)  boundaries. 

Based  on  the  above  approach,  a  theory  was  formulated  for  finding 
the  pressure  distribution  on  any  airfoil  at  angle  of  attack  in  a  nonuniformly 


7 


sheared  flow  with  a  piecewise -linear  velocity  profile.  The  computer 
program  has  been  written  for  a  two-dimensional  velocity  profile  like 
that  of  Figure  1(b).  This  velocity  profile  has  the  advantage  of  being 
the  simplest  nonuniformly  sheared  flow  to  which  the  theory  can  be 
applied  ;  this  is  a  worthwhile  consideration  for  the  initial  calcu¬ 
lations.  Details  of  the  theory  and  its  implementation  on  the  CAL  IBM 
7044  digital  computer  are  presented  in  the  following  sections  of  this 
report. 


It  is  difficult  to  assess  with  any  exactness  the  implications  of  the 
restriction  to  piecewise-linear  velocity  distributions.  In  effect,  a  free- 
stream  vorticity  distribution  like  that  of  the  solid  line  in  the  following 
sketch  has  been  replaced  with  one  like  that  of  the  dashed  line. 


VORTICITY 


Perhaps  the  best  way  to  consider  this  aspect  is  that  the  theory 
represents  a  limiting  case  of  infinite  shear  gradient  and  that  this  is  the 
only  restriction  on  the  model.  It  might  be  that  the  ratio  of  the  distance  £_ 
to  airfoil  thickness  as  well  as  the  ratio  of  X  to  the  displacement  of  the 
airfoil  chord  away  from  the  discontinuity  must  be  small  for  the  approxi¬ 
mation  to  be  valid. 
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DESCRIPTION  OF  COMPUTER  PROGRAM 


Details  of  the  formulation  and  implementation  of  the  computer 
program  are  presented  in  the  appendix.  A  brief  outline  of  the  basic 
ideas  incorporated  and  certain  numerical  techniques  of  importance  to 
the  program  are  discussed  in  the  following  paragraphs. 

Figure  2  is  a  schematic  drawing  of  the  two-dimensional  flow 
model  assumed  for  the  computer  program.  A  Cartesian  coordinate 
system  is  used.  The  *  -axis  is  parallel  to  the  wind-tunnel  walls  and 
positive  in  the  direction  of  the  flow.  The  origin  is  centered  vertically 
between  the  walls, and  %  =  0  is  located  at  the  airfoil  midchord.  The 
wind-tunnel  walls  and  the  airfoil  surface  are  represented  by  bound 
vortex  sheets  of  varying  strength.  The  free-stream  velocity  variation 
in  the  y -direction  is  piecewise  linear  and  symmetric  about  y  -  0  (point 
(a)  in  Figure  2).  On  either  side  of  the  streamline  which  passes  through 
point  (a)  (denoted  henceforth  as  the  "dividing  streamline"),  the  distributed 
vorticity  is  constant.  If  the  assumed  free-stream  velocity  profile,  U(y)  , 
is  symmetric  {  Uf  =  U2  in  Figure  2),  as  is  the  case  in  the  work  presented 
here,  the  vorticity  on  either  side  of  the  dividing  streamline  is  equal  in 
magnitude  but  of  opposite  sign.  If  the  airfoil  is  of  nonsymmetric  profile, 
or  is  nonsymmetrically  disposed  in  the  wind  tunnel  (for  example,  at 
zero  angle  of  attack  but  not  centered  between  the  walls,  or  at  angle  of 
attack),  the  resulting  disturbance  flow  will  distort  the  dividing  stream¬ 
line  in  some  fashion.  This  distorted  shape  is  denoted  by  gU)  in  Figure  2. 
Equation  (1)  states  that  the  vorticity  in  a  rotational  flow  is  constant 
along  streamlines;  therefore,  the  dividing  streamline  still  separates 
regions  of  constant  vorticity.  If  the  shape  of  the  dividing  streamline 
q(z)  were  known,  then  it  would  be  possible  to  compute  the  variation  in 
strength  of  the  bound  vortex  sheets  representing  the  wind-tunnel  walls 
and  airfoil  surface.  Given  this  information,  plus  information  locating 
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the  position  of  the  stagnation  streamline  which  intersects  the  airfoil, 
the  pressure  distribution  on  the  airfoil  surface  can  be  computed. 

An  important  part  of  the  computer  program  is  an  iterative 
technique  for  determining  the  function  g(x)  numerically.  Mass -flow 
continuity  is  the  basis  for  this  iterative  technique  ir.  the  developed 
computer  program.  If  an  initial  shape  y0(x>)  is  assumed,  the  strengths 
of  the  singularity  distributions  (vortex  sheets  representing  the  wind- 
tunnel  walls  and  the  airfoil  surface)  can  be  determined.  Once  these 
distributions  are  known,  the  first  iterated  streamline  shape,  ,  is 

determined  by  first  computing,  at  each  value  of  x  considered,  the 
vertical  point  (  qf(x)  ,  say)  in  the  flow  where  the  mass  flow  is  equal  to 
that  on  either  side  of  the  dividing  streamline  in  the  undisturbed  free 
stream,  and  then  by  taking  the  average  of  Q0(x)  and  q*(x)  .  Thus 

1t(»)  -  + ?/<*)]• 

Given  g,(x)  ,  again  the  singularity  strength  distributions  for  the  bound 
vortex  sheets,  representing  the  wind-tunnel  walls  and  the  airfoil,can 
be  obtained, from  which  the  second  iterated  streamline  shape,  qg(x>)  , 
can  be  determined.  If  the  iterative  process  converges,  Um  gt(x)  -  g(x). 
Hopefully,  the  rate  of  convergence  is  such  that  this  limit  is  very  nearly 
attained  after  only  a  few  iterations. 

The  singularity  strength  distributions  required  to  satisfy  the 
boundary  conditions  at  the  wind-tunnel  walls  and  the  airfoil  surface  are 
those  which  correspond  to  the  g(x)  evaluated  in  the  final  iteration.  The 
flow  field  is  thus  completely  determined.  Velocities  on  the  airfoil 
surface  can  then  be  converted  to  static  pressures  by  means  of  the 
Bernoulli  equation  for  rotational  flow, 

■Yfi(u*+v‘)+p  =  3 />  [uf  !/)]*  + ft,  .  (2) 


10 


applied  along  the  airfoil  stagnation  streamline.  The  proper  value  of  U(y) 
in  the  free  stream  to  be  used  in  Equation  (2)  is  determined  by  computing 
the  mass  flow  between  the  wind-tunnel  wall  and  the  airfoil.  The  value 
of  U(y)  which  corresponds  to  this  mass  flow  in  the  undisturbed  free 
stream  is  the  correct  value. 


DISCUSSION  OF  COMPUTER  PROGRAM 


The  strengths  of  the  singularity  distributions  (vortex  sheets) 
representing  the  wind-tunnel  walls  and  airfoil  surface  must  be  evaluated 
in  order  to  be  able  to  calculate  airfoil  pressure  distributions.  This  is 
done  in  the  computer  program  by  breaking  up  the  vortex  sheets  into 
small  segments.  For  each  of  these  segments,  the  singularity  strength 
is  assumed  constant.  The  integrals  over  each  segment  of  the  vortex 
sheets  for  use  in  the  equations  for  velocity  components  (Equations  (A -2), 
(A-3),  (A-6),  and  (A-7)  of  the  appendix  )  can  then  be  evaluated  in  closed 
form.  Hence,  the  velocity  components  are  expressible  as  finite  series 
in  terms  of  the  vortex  sheet  singularity  strengths.  The  vortex  sheets 
representing  the  wind-tunnel  walls  are  assumed  to  have  known  constant 
strengths  upstream  and  downstream  of  the  airfoil.  The  boundary 
conditions  for  no  normal  flow  across  the  wind-tunnel  walls  and  airfoil 
surface  (Equations  (A-8)  and  (A-9)  of  the  appendix)  are  next  applied  at 
the  center  of  each  of  the  vortex  sheet  segments.  If  these  boundary 
conditions  were  the  only  conditions  to  be  satisfied,  a  determinate  set 
of  linear  simultaneous  algebraic  equations  would  result.  However, 
in  addition,  the  Kutta-Joukowski  condition  at  the  airfoil  trailing  edge 
must  be  satisfied.  This  constraint  is  implemented  in  the  computer 
program  by  requiring  that  the  sum  of  the  singularities  representing  the 
upper  and  lower  airfoil  surfaces  vanish  in  the  limit  as  the  trailing  edge 
is  approached.  This  results  in  an  additional  equation,  and  the  problem 
is  now  overdetermined.  The  problem  can  be  made  determinate  by 
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relaxing  the  flow-tangency  requirement  at  one  of  the  vortex  sheet 
segments.  However,  additional  considerations  arose  in  the  solution  of 
the  set  of  simultaneous  equations  which  obviated  this  necessity. 

Clearly,  if  the  problem  is  properly  posed  in  a  mathematical  sense,  the 
equations  which  result  should  be  determinate.  The  overdeterminacy  of 
the  above  problem  with  the  added  Kutta- Joukowski  condition  most  likely 
results  from  the  method  of  numerical  analysis  used  to  obtain  solutions. 

For  example,  if  it  had  been  assumed  that  the  singularity  strength  of 
each  segment  of  vortex  sheet  varied  linearly,  there  would  have  been  as 
many  equations  as  unknowns.  The  choice  of  constant  singularity 
strength  for  each  segment,  rather  than  a  linear  variation,  was  dictated 
primarily  by  the  necessity  to  limit  the  number  of  linear  equations  which 
had  to  be  solved. 

Effectively,  the  process  described  above  for  determining  the 
vortex  strength  distributions  is  the  numerical  solution  of  an  integral 
equation.  One  check  on  the  computer  program  was  the  computation  of 
the  pressure  distribution  on  the  airfoil  in  a  uniform  flow.  As  the  air¬ 
foil  involved  is  a  Joukowski  17-percent-thick  symmetric  airfoil,  analytic 
solutions  for  uniform  inviscid,  two-dimensional  flow  are  available. 

Initial  numerical  calculations  for  the  airfoil  at  an  angle  of  attack 
of  10  degrees  resulted  in  extremely  irregular  variation  with  %  of 
computed  values  of  Cp  on  the  airfoil  surface,  as  shown  in  Figure  3. 
Reference  10  discusses  similar  behavior  exhibited  by  numerically 
computed  solutions  to  integral  equations  of  the  first  kind  and  presents  a 
technique  for  obtaining  smoothly  varying  solutions  in  an  optimum  sense, 

where  it  is  known  beforehand  that  the  solution  should  be  smooth.  Reference 
11  also  considers  the  same  problem,  and  extends  the  smoothing  technique 

to  allow  the  treatment  of  overdetermined  sets  of  simultaneous  algebraic 
equations.  As  previously  noted,  when  the  Kutta- Joukowski  condition 

is  included  in  the  present  computer  program,  an  overdetermined  set 
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of  equations  is  obtained.  Hence,  the  modified  smoothing  technique  of 
Reference  11  was  incorporated  into  the  solution  for  the  vortex  strength 
distributions  of  the  computer  program;  details  are  presented  in  the 
appendix.  Figure  4  shows  computed  values  of  Cp  on  the  airfoil  surface 
at  a  -  10  degrees,  which  resulted  when  the  smoothing  technique  was 
used  in  the  computer  program.  It  is  evident  that  excellent  agreement 
with  the  exact  solution  to  the  uniform,  inviscid  flow  problem  was  obtained, 
with  the  exception  of  a  3-percent  overshoot  at  the  maximum  negative  com¬ 
puted  value  of  Cp  (suction  peak).  Apparently,  the  smoothing  process  works 
quite  well  in  the  present  instance. 

In  contrast  to  the  uniform  flow  calculation  noted  above,  the 
nonuniformly  sheared  flow  calculations  required  an  iterative  process  to 
obtain  solutions.  Although  there  was  confidence  (based  on  the  uniform 
impinging  jet  solutions  obtained  previously  and  presented  in  Reference  9) 
that  a  convergent  iterative  technique  could  be  derived,  the  development 
of  one  which  converged  much  more  rapidly  than  that  for  the  impinging 
jet  constituted  a  major  objective  of  the  present  work.  In  that  previous 
case, as  many  as  80  iterations  were  required  to  obtain  a  satisfactory  solu¬ 
tion.  Since  the  present  computer  program  required  about  12  minutes  of 
machine  time  per  iteration,  80  iterations  would  represent  an  unacceptable 
16  hours  of  computer  time.  The  technique  finally  employed  was  successful 
in  reducing  the  number  of  iterations  required;  at  most,  five  to  six  were 
sufficient  to  insure  an  acceptable  solution  for  the  configurations  run  thus 
far.  This  represents  an  hour  to  an  hour  and  fifteen  minutes  of  IBM  7044 
computer  time. 

It  should  be  noted  that  the  use  of  mass  flow  continuity  in  the 
iterative  process  requires  a  known  reference  streamline  of  the  flow.  In 
the  present  case,  a  wind-tunnel  wall  provides  this  reference  streamline. 
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EXPERIMENTAL  PROGRAM 


GENERAL 

The  experimental  program  was  designed  to  provide  airfoil 
pressure  distribution  data  from  which  the  accuracy  of  the  theoretical 
results  could  be  determined.  At  the  same  time,  additional  experi¬ 
mental  evidence  regarding  airfoil  maximum  lift  behavior  in  sheared 
flows  would  be  obtained. 

The  wind-tunnel  airfoil  tests  were  of  three  types: 

1.  Uniform  flow,  low-turbulence  tests 

2.  Uniform  flow,  high-turbulence  tests 

3.  Sheared  flow  tests 

The  velocity  profile  generated  for  the  wind-tunnel  tests  of  the  airfoil 
in  sheared  flow  was  similar  to  the  particular  profile  considered  in  the 
theoretical  treatment  of  the  problem  during  this  program. 

The  uniform-flow  ests,  with  and  without  high  turbulence,  were 
designed  as  an  experimental  control.  The  shear  screen  used  to  generate 
the  sheared  flow  also  generated  a  high  turbulence  level  in  the  flow.  It 
was  possible  that  this  high  turbulence  level  would  in  itself  have  had  some 
marked  effect  on  airfoil  separation  characteristics.  The  low  and  high 
turbulence  tests  in  uniform  flow  were  to  resolve  this  possibility.  In  the 
event  that  the  effects  of  turbulence  in  uniform  flow  were  small,  then  the 
high-turbulence,  uniform-flow  airfoil  pressure  distribution  data  would 
provide  a  reference  insofar  as  a  comparison  of  theory  and  experiment 
is  concerned.  A  comparison  of  theoretical  and  experimental  pressure 
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distributions  for  the  sheared-flow  case  could  not  be  expected  to  be  in 
any  better  agreement  than  the  pressure  distributions  from  uniform  flow 
theory  and  experiment.  The  boundary  layer  is  known  to  alter,  somewhat, 
the  experimental  pressure  distribution  when  compared  to  theory,  even 
at  angles  of  attack  which  are  significantly  below  the  stall  angle. 

EXPERIMENTAL  APPARATUS 

The  experiments  were  made  in  the  subsonic  leg  of  the  CAL/AF 
One -Foot  High-Speed  Wind  Tunnel.  This  leg  of  the  wind  tunnel  has  a 
test  section  with  a  cross  section  of  17  inches  by  24  inches  and  is  operated 
as  a  closed-throat  nonreturn  tunnel.  With  no  shear  screen,  the  test 
section  stagnation  pressure  is  one  atmosphere,  and  the  practical  speed 
range  in  the  clear  test  section  is  from  0  to  approximately  120  feet  per 
second. 


The  two-dimensional  airfoil  used  in  this  research  has  a  symmetric 
Joukowski  profile  with  a  thickness-chord  ratio  of  17  percent  and  a  chord  of 
6  inches.  A  number  of  static  pressure  orifices  were  distributed  chordwise 
on  the  top  and  bottom  surfaces  along  a  section  Located  near  the  center  of 
the  wing  span.  The  wing  was  mounted  in  the  wind-tunnel  test  section  by  a 
supporting  structure  which  permitted  independent  change  in  geometric  angle 
of  attack  and  vertical  location  of  the  wing  in  the  wind  tunnel  by  means  of 
manual  adjustment.  Figure  5  shows  the  wing  and  wing -support  system 
mounted  in  the  wind-tunnel  test  section. 

Instrumentation  was  relatively  simple.  Flow  calibrations  in 
the  wind-tunnel  test  section  were  performed  with  a  conventional  3/16- 
inch-diameter  Pitot-static  probe  which  included  static -pressure  taps 
from  which  flow  angularity  data  were  obtained.  Static  and  dynamic 
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pressures  were  measured  with  manometers.  Thirty-two  static-pressure 
taps  distributed  over  the  top  and  bottom  surfaces  of  the  wing  were  con¬ 
nected  to  an  inclined  manometer  bank,  and  all  manometer  readings  were 
recorded  photographically.  Free-stream  static  pressure  was  obtained 
from  a  pressure  tap  on  the  ceiling  of  the  test  section.  This  tap  was  approxi¬ 
mately  2-1/2  chord  lengths  upstream  of  the  model  midchord. 

The  sheared  flow  in  the  wind-tunnel  test  section  was  produced 
by  a  screen  placed  slightly  more  than  3  feet  (between  6  and  7  wing  chords) 
upstream  of  the  wing.  This  screen,  shown  in  Figure  6,  consisted  of  a 
number  of  circular  rods  spanning  the  wind-tunnel  section  horizontally 
and  secured  by  a  frame  which  was  clamped  between  two  sections  of  the 
wind-tunnel  circuit.  The  spacing  between  rods  and  the  rod  diameters 
were  varied  so  as  to  introduce  variable  losses  across  the  flow.  By 
proper  spacing  and  rod  size  variation,  the  vertical  distribution  of  losses 
in  the  flow  at  the  screen  can  be  such  that  the  desired  sheared-flow  velocity 
contour  can  be  obtained  in  the  test  section. 

The  theory  presented  in  Appendix  II  of  Reference  5  was  used  in 
the  design  of  the  screen.  Application  of  this  theory  is  somewhat 
laborious,  and  the  results  are  dependent  on  the  accuracy  of  the  empirical 
relationship  between  local  screen  solidity  and  loss  coefficient. 

Apparently  the  empirical  curve  presented  in  Reference  5  is  not  reliable 
at  high  values  of  solidity,  because  the  use  of  this  curve  provided  a  velocity 
distribution  which  was  considerably  different  from  the  design  velocity 
distribution  near  the  centerline  of  the  test  section.  The  final 
configuration  used  for  the  shear  screen  was  tailored  by  a  trial  and 
error  process  of  small  modifications  to  the  theoretical  design  until  it 
was  judged  that  further  improvement  could  not  be  obtained  without 
undue  expenditure  of  time. 
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WIND-TUNNEL  TESTS 


Airfoil  pressure  distributions  were  measured  in  the  wind  tunnel 
with  uniform  flow  and  with  a  two-dimensional,  nonuniformly  sheared 
flow  in  the  test  section.  Uniform-flow  airfoil  data  were  obtained  with 
ambient  turbulence  levels  (clear  tunnel)  ana  with  high  turbulence  levels 
generated  by  screens  mounted  upstream  of  the  test  section.  These 
screens  were  similar  to  the  shear  screen  except  that  the  horizontal 
rods  used  were  uniformly  spaced.  Two  such  turbulence  screens  were 
used  in  separate  tests,  one  screen  with  rod  diameters  and  spacing  corres¬ 
ponding  to  that  portion  of  the  shear  screen  with  low  solidity  (relatively 
large  open  area  ratio)  and  one  screen  corresponding  to  the  high  solidity 
portion  of  the  shear  screen. 

During  the  uniform -flow  tests,  data  were  obtained  with  the  wing 

model  spanning  the  test  section  horizontally  midway  between  the  top 

and  bottom  walls.  Pressure  data  were  recorded  within  a  range  of  angles 

of  attack  from  near  zero  lift  to  beyond  maximum  lift.  Two  wind  speeds 

55 

were  used  to  produce  Reynolds  numbers  of  1.8  x  10  and  3.  7  x  10  ,  based 
on  the  airfoil  chord  of  6  inches. 

Wind-tunnel  flow  calibrations  were  performed  for  both  the  clear - 
tunnel  uniform-flow  condition  and  with  the  turbulence  screens  in  place. 

In  the  clear  tunnel,  variation  of  dynamic  pressure  in  the  free  stream 
at  the  model  location  was  less  than  1  percent  over  the  model  chord, 
and  variation  in  flow  angularity  in  the  vertical  plane  was  less  than  0.  4 
degree  at  the  highest  speed  tested  and  about  0,  1  degree  at  the  lowest 
speed.  With  the  low-solidity,  uniform  screen  installed,  vertical  traverses 
of  the  flow  at  the  mode]  midchord  reference  station  in  the  test  section 
resulted  in  maximum  variation  of  dynamic  pressure  of  ±  2.  8  percent 
and  maximum  flow  angularity  variation  in  the  vertical  plane  of  ±  0.  25 
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degree.  The  flow  with  the  high-solidity,  uniform  screen  installed  was 
considerably  less  uniform,  with  maximum  variations  in  free-stream 
dynamic  pressure  of  as  much  as  ±  19  percent.  A  portion  of  the  flow 
from  the  test  section  centerline  to  4  inches  below  was  more  nearly 
uniform,  and  for  wind-tunnel  tests  with  this  turbulence  screen,  the 
model  was  centered  in  this  region. 

Sheared-flow  test  data  were  obtained  with  the  wing  model 
situated  at  varying  heights  above  and  below  the  centerline  of  the 
sheared  flow.  At  each  position  of  the  wing  in  the  test  section, 
pressure  distributions  were  measured  within  a  range  of  angles  of 
attack  from  near  zero  lift  to  beyond  maximum  lift.  The  nominal 
Reynolds  number,  based  on  the  airfoil  chord  and  velocity  at  the  posi¬ 
tion  of  the  midchord,  was  different  for  each  height  because  of  the 

variation  of  velocity  with  height.  The  range  of  Reynolds  numbers 

5  5 

was  from  1.7x10  at  the  tunnel  centerline  to  2,  5  x  10  at  3  inches 

above  and  below  the  centerline.  From  previous  experience  with  the 

wing  model  and  wind  tunnel,  it  was  known  that  the  boundary  layers  on 

the  wind-tunnel  side  walls  tended  to  interact  with  the  wing  model  boundary 

layer,  particularly  near  stall,  so  that  stall  occurred  in  a  three*  iimen- 

sional  manner.  Wind-tunnel  wall  suction  near  the  intersection  of  the 

wing  and  the  wall  was  used  to  minimize  this  interaction.  Optimum 

values  of  the  wall  suction  were  obtained  by  means  of  lampblack  and 

kerosene  flow  visualization  studies,  which  preceded  pressure  tests 

for  each  configuration  change. 

Velocity  profiles  generated  by  the  shear  screen  were  obtained 
through  a  series  of  flow  calibrations  in  the  test  section  in  a  vertical 
plane  at  the  model  reference  axis  (midchord  of  the  wing),  using  the 
3/16-inch-diameter  Pitot-static  probe.  Figure  7  presents  velocity 
distributions  for  the  final  shear  screen  configuration  used  during  pres¬ 
sure-data  tests.  Wing  test  data  were  obtained  with  the  pressure-tapped 
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wing  model  segment  centered  0.  75  inch  to  one  side  (left  looking  down¬ 
stream)  of  the  wind-tunnel  centerline.  It  can  be  seen  that  the  velocity 
distribution  approached  the  design  distribution  with  reasonable  accuracy 
but  that  some  irregularities  were  present  in  the  experimental  data  in 
the  vicinity  of  5  inches  above  and  below  the  tunnel  centerline.  These 
irregularities  are  not  believed  to  be  of  importance,  because  they  were 
not  large  and  because  the  test  region  extended  only  to  3  inches  above  and 
below  the  centerline.  Of  perhaps  more  concern  is  the  rounded  shape  of 
the  experimental  velocity  profile  near  the  centerline.  This  rounded  shape 
is  not  in  agreement  with  the  mathematical  model  and  will  be  shown  later 
to  cause  uncertainties  in  the  data  reduction  technique  used  to  define  the 
stagnation  streamline  which  intersects  the  wing.  It  is  noted,  however, 
that  such  profiles  are  more  likely  to  be  encountered  in  practice. 

Figure  8  shows  the  flow  angularity  at  the  plane  of  the  model 

midchord  as  determined  with  the  static -pressure  angularity  taps  on  the 

3/ 16 -inch  probe.  The  probe  was  calibrated  for  these  measurements  using 

standard  procedures  as  described,  for  example,  in  Reference  12.  With 

the  probe  tip  fixed  at  some  location  in  the  flow,  the  probe  is  rotated  through 

an  arc  in  the  vertical  plane,  say  ±  10  degrees  from  horizontal.  The  probe 

is  then  inverted  and,  with  the  probe  tip  at  precisely  the  same  location  in 

the  flow,  the  procedure  is  repeated.  The  output  signals  are  plotted  against 

angular  position  for  each  case.  The  slopes  and  the  crossing  point  of  the 

two  sets  of  data  provide  the  sensitivity  of  the  probe  for  measuring  flow 

angularity  and  a  correction  for  the  probe's  dissymmetry  which  are  inde- 

* 

pendent  of  the  characteristics  of  the  flow  itself.  The  calibration  was 
checked  at  several  vertical  locations  in  the  sheared  flow  even  though  it 
was  expected  to  be  unaffected  by  the  flow.  It  was  found,  however,  that 


*  The  data  reduction  included  a  small  correction  for  the  theoretical 
effect  of  shear  on  the  measurements. 
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the  calibration  was  valid  only  in  regions  far  removed  from  the  centerline 
of  the  sheared  flow  (  |y|  >  3  inches).  In  the  vicinity  of  the  center  of  the 
sheared  flow,  the  sensitivity  of  the  probe  remained  constant  (as  in  the 
outer  region)  but  the  required  correction  for  probe  dissymmetry  appeared 
to  vary.  (This  was  possibly  due  to  the  effects  of  curvature  of  the  velocity 
profile  in  this  central  region  of  the  flow. )  Measurement  of  the  flow  angu¬ 
larity  with  sufficient  precision  required,  therefore,  that  a  technique 
similar  to  the  probe-calibration  procedure  be  used  near  the  centerline 
of  the  sheared  flow.  At  each  vertical  location  within  3  inches  from  the 
centerline  at  which  flow  angularity  was  measured,  the  probe  was  rotated 
through  an  arc  in  the  vertical  plane  with  the  probe  both  right-side  up  and 
inverted.  The  crossing  point  of  the  two  sets  of  data  provided  the  measure 
of  the  flow  angularity  independent  of  the  apparent  correction  for  probe 
dissymmetry.  For  locations  outside  of  3  inches  from  the  centerline,  the 
flow  angularity  could  be  determined  by  the  substantially  simpler  method 
of  using  the  probe  calibration  and  the  yaw-tap  measurements  made  while 
the  probe  was  traversed  vertically  with  a  fixed  geometric  an^le  for  its 
axis.  It  is  for  this  reason  that  Figure  8  shows  a  much  higher  density  of 
data  points  in  the  outer  region  than  in  the  central  portion  of  the  sheared 
flow. 

A  mean  curve  of  flow  angularity  versus  vertical  location  which 
was  fitted  to  the  data,  as  shown  in  Figure  8,  has  a  variation  of  only  from 
-0.  30  degree  to  +  0.  35  degree.  Accuracy  of  this  mean  curve,  based  on 
the  repeatability  of  the  measurements,  is  judged  to  be  within  ±0.  25 
degree.  The  mean  curve  was  used  in  calculating  the  effective  angle 
of  attack  of  the  wing  for  each  vertical  location. 
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DISCUSSION  OF  EXPERIMENTAL  RESULTS 


Although  a  very  large  number  of  pressure  distributions  (over 
200)were  measured  during  the  course  of  this  program,  only  distri¬ 
butions  which  serve  to  illustrate  the  points  under  discussion  will  be 
presented.  The  section  lift  coefficient,  ,  provides  a  convenient  method 
for  identifying  regions  where  unusual  phenomena  may  be  occurring. 

Hence,  in  this  report,  curves  of  c ^  versus  angle  of  attack,  oc  ,  will  be 
presented  first.  This  coefficient,  along  with  section  drag  and  moment 
coefficients,  was  obtained  by  numerical  integration  of  the  pressure 
coefficient  distribution  on  the  wing  section. 


UNIFORM  FLOW  TESTS 


Figure  9  shows  Cj  versus  a  for  the  Joukowski  airfoil  under  a 
variety  of  uniform-flow  test  conditions.  These  test  conditions  include 
tests  at  different  free-stream  velocities  and  tests  with  and  without  the 
low-solidity  turbulence  screen.  The  classical  theoretical  curve  for  a 
17-percent-thick  symmetrical  Joukowski  airfoil  and  a  curve  obtained 
from  the  strain-gauge  balance  measurements  of  Reference  5  on  the 
same  airfoil  are  shown  for  comparison.  Data  points  obtained  after 
the  airfoil  had  stalled  suddenly  and  completely  are  noted.  It  can  be 
seen  that  the  experimental  data,  including  that  of  Reference  5,  are 
in  essential  agreement,  with  the  possible  exception  of  the  clear  tunnel 
data  point  at  00  =  1 1 8  f,  p.  s.  and  a  =  7,  5  degrees. 

The  maximum  value  of  Cj  attained  and  the  values  of  Cj  at  angles 
larger  than  that  for  the  maximum  appear  to  be  more  dependent  on  the 
free-stream  velocity  than  on  whether  the  flow  was  passed  through  a 
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screen  before  encountering  the  airfoil.  In  general,  the  higher  velocity 
free  stream  had  a  lower  value  of  Ci  for  a  given  oc.  beyond  that  for 
maximum  Cj  .  The  overall  difference,  however,  is  small.  At  angles 
below  that  for  maximum  Cj  ,  the  trend  appears  to  be  reversed.  All 
the  experimental  data  fall  below  the  theoretical  curve  in  a  fashion 
similar  to  that  expected  from  experience.  Data  obtained  with  the  high- 
solidity  turbulence  screen  in  place  are  not  shown  because  of  difficulties 
encountered  in  measuring  the  angularity  of  the  flow.  The  maximum 

value  of  Cj  which  was  obtained  with  this  screen  was  1.  25,  which  is  in 
agreement  with  the  data  of  Figure  9. 

Pressure  distributions  for  angles  of  0,  7.  5,  and  10.0  degrees  are 
shown  in  Figures  10,  11,  and  12  and  are  compared  with  the  inviscid  poten¬ 
tial  flow  theory.  The  agreement  between  theory  and  experiment  is  satis¬ 
factory  at  os  =  0  degrees,  but  it  begins  to  deteriorate  near  the  upper  surface 
suction  peak  and  near  the  trailing  edge  as  (X  is  increased.  The  presence 
of  the  turbulence  screen  in  the  tunnel  lowers  the  peak  suction  attainable 
on  the  upper  surface;  but  it  appears  to  delay  the  onset  of  separation  from 
the  upper  surface  near  the  trailing  edge.  This  separation  is  indicated 
by  a  leveling  of  the  value  of  Cp  near  the  trailing  edge  and  is  mildly 
discernable  in  Figure  12  near  t/c  »  0.  85  for  the  clear  tunnel  pressure 
distribution. 

There  are  apparent  discontinuities  in  the  experimental  data  for 
the  upper  surface  aft  of  the  suction  peak;  for  example;  between  */c  =  0.  1  5 
and  0.  20  on  Figures  11  and  12  for  the  clear  tunnel  data.  The  flow 
visualization  studies  showed  that  these  irregularities  correspond  to  the 
turbulent  reattachment  of  a  leading- edge  separation  bubble.  Such 
bubbles  occurred  even  at  cc,  -  0  degrees.  They  decreased  in  length  as  \ 

the  angle  of  attack  was  increased  or  the  flow  velocity  was  increased. 

The  addition  of  the  turbulence  screen  to  the  flow  did  not  give  a 
measurably  consistent  variation  to  the  bubble  length.  According  to 
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Reference  10,  which  is  an  excellent  survey  of  such  bubbles,  these 
laminar  separation  bubbles  would  generally  be  classed  as  short  bubbles, 
and  are  evidently  relatively  thin. 

Pressure  distributions  obtained  on  the  airfoil  near  maximum 
lift  and  after  complete  stall  are  shown  in  Figures  13,  14,  and  15  for 
66  f.  p.  s.  and  with  the  low-solidity  screen  installed.  This  free-stream 
velocity  with  the  turbulence  screen  installed  corresponds  most  closely 
to  the  conditions  obtained  when  the  shear  screen  was  used.  Figures  13 
and  14  clearly  show  the  progressive  forward  movement  of  the  separation 
point  on  the  upper  surface  and  also  show  the  forward  movement  of  the 
trailing  edge  of  the  bubble.  At  a  =  13.0  degrees,  which  corresponds  to 
maximum  Cj  for  this  case,  separation  occurs  near  i/c  -  0.  65.  Between 
<2  =  13.  0  and  16.  5  degrees, the  suction  peak  remains  constant  (Cp-  -9.Z) , 
but  separation  occurs  near  z/c  =  0.  28  for  the  larger  angle  of  attack. 

Full  separation  or  stall  occurred  at  &  -  16.  7  degrees,  as  is  evident  in 
Figure  15. 

Figures  16  and  17  show  respectively  the  values  of  cm  and  cd 
versus  obtained  with  the  low-solidity  screen  in  place.  The  data  shown 
in  Figures  16  and  17  must  be  regarded  with  some  reserve,  since  they  are 
obtained  from  the  small  difference  between  two  large  force  components. 
The  consistency  of  the  data,  however,  suggests  that  the  accuracy  was  at 
least  as  good  as  that  obtained  with  the  strain-gauge  balance  of  Reference  5. 
The  drag  coefficient  data  of  Figure  17  were  derived  from  pressure  distri¬ 
butions;  hence,  skin  friction  is  not  included  in  cd. 

In  general,  the  uniform  flow  data  provided  the  expected  results, 
although  the  existence  of  the  leading -edge  bubble  at  low  angles  of  attack 
was  not  expected.  The  addition  of  the  turbulence  screen  did  influence 
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the  data  very  slightly,  but  not  to  the  extent  necessary  to  provide  an 
explanation  of  results  observed  with  the  shear  screen  in  place.  These 
latter  results  are  presented  in  the  next  section.  Since  there  was  a 
slight  dependence  on  whether  or  not  the  turbulence  screen  was  in  place, 
the  uniform  flow  data  used  for  a  control  should  be  that  obtained  with  the 
turbulence  screen  in  place.  The  appropriate  free-stream  velocity  is 
U0  =  66  f.  p.  s.  ,  since  this  is  approximately  midway  between  the  minimum 
and  maximum  velocities  encountered  with  the  shear  screen.  From  the 
comparison  of  the  measured  and  the  theoretical  pressure  distributions 
for  uniform  flow,  an  increasing  divergence  between  theory  and  experiment 
for  the  sheared  flow  can  be  expected  as  angle  of  attack  is  increased,  with 
the  theory  predicting  higher  suction  peaks  on  the  upper  surface  than  are 
evidenced  in  the  corresponding  experimental  data. 


SHEARED  FLOW  TESTS 


The  research  reported  here  had  as  its  primary  objective  the 
completion  of  the  development  of  the  computer  program  to  compute 
airfoil  characteristics  in  the  nonuniformly  sheared  flow  of  Figure  1(b)  and 
the  validation  of  the  theory  by  means  of  comparison  with  experiment. 

An  added  objective  was  to  investigate  experimentally  in  some  detail  the 
aerodynamic  behavior  and,  particularly,  the  stalling  characteristics  of 
an  airfoil  in  a  nonuniformly  sheared  flow. 

The  significant  feature  of  the  nonuniformly  sheared  flow  velocity 
profile  used  in  the  present  work  is  that  there  are  no  free  jet  boundaries. 

It  would  seem  that  if  an  airfoil  intersects  a  nonuniformly  sheared  flow 
which  is  wholly  contained  within  a  jet  with  free  boundaries,  as  in  Figure 
1(c),  the  deformation  of  these  free  boundaries  might  have  significant 
effects  on  the  aerodynamic  characteristics.  Hence,  the  choice  of  a 
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nonuniformly  sheared  flow  without  free  jet  boundaries,  besides 
simplifying  the  mathematical  problem, allowed  experimental  isolation  of 
the  effect  of  nonuniform  shear  on  airfoil  aerodynamics  from  the  effect 
of  free  jet  boundaries. 

A  major  problem  in  discussing  the  aerodynamics  of  airfoils  in 
sheared  ilows  is  the  question  of  a  proper  reference  dynamic  pressure. 
The  midchord  free-stream  dynamic  pressure  is  convenient  for  data 
reduction.  From  a  practical  viewpoint,  the  lift  generated  is  of  most 
interest;  hence,  some  constant  value  of  dynamic  pressure  should  be  used. 
However,  in  terms  of  boundary  layer  characteristics,  including 
separation,  the  actual  surface  static  pressures  and  associated  pressure 
gradients  are  of  primary  importance;  hence,  the  stagnation  streamline 
free-stream  dynamic  pressure  is  the  best  reference  dynamic  pressure. 

In  fact,  it  is  this  multiplicity  of  possible  reference  dynamic  pressures  in 
a  sheared  flow  which  tends  to  cause  confusion,  and  it  possibly  leads  to 
misleading  conclusions  regarding  experimental  aerodynamic  data. 

The  free-stream  dynamic  pressure  of  the  stagnation  streamline 
is  equal  to  the  difference  between  the  static  pressure  at  the  stagnation 
point  on  the  airfoil  and  the  free-stream  static  pressure.  J.f  taps  are 
spaced  sufficiently  close  together  on  the  lower  leading  edge  surface  of 
the  airfoil,  it  is  possible  to  obtain  the  stagnation  streamline  free-stream 
dynamic  pressure  from  the  airfoil  pressure  distributions.  Bernoulli's 
equation  along  the  stagnation  streamline  is 

p  -h  ~p(uZt  vz)  -  p0  +  J  p(^s)2 

where  pa  is  the  free-stream  static  pressure  and  l^is  the  velocity  in  the 
sheared  free  stream  which  corresponds  to  the  stagnation  streamline. 
Using  this  relation,  the  pressure  coefficient,  ( Cp)f ,  based  on  the  free- 
stream  dynamic  pressure  of  the  stagnation  streamline  becomes 
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P-Po  U*+V2 

'c')s  *ir<usf  = 

Note  that  since  the  maximum  positive  value  of  (C/t)i  is  1,  the  appro¬ 
priate  free-stream  dynamic  pressure  of  the  stagnation  streamline  is 
simply  the  maximum  positive  value  of  (p-p0)  which  is  reached  on  the 
airfoil.  In  addition,  the  remaining  two  reference  dynamic  pressures 
which  were  mentioned  above  provide  pressure  coefficients  Cp  (based  on 
an  average  velocity,  U)  ,  and  (Cp)R  (based  on  midchord  velocity,  UR  ), 
which  are  defined  by 

-  _  p' Po  ( ) 1  uZ+  w  * 

Cp  "  =  \0  >  (0)z 

,r\  P‘Po  /i£ )* 

[  p]^{p(u*)2  \uj  '  (UP2 

The  section  force  and  moment  coefficients  can  be  referred  to  any  of  the 
reference  dynamic  pressures  to  obtain,  for  example,  ,  (cpR,  or  (cx)s  . 

Figure  18  shows  the  variations  of  average  lift  coefficient,  c/  , 
with  angle  of  attack  which  were  obtained  in  the  present  investigation. 

In  this  case,  Cj  is  referred  to  a  dynamic  pressure  derived  from  the 
average  free-stream  velocity  in  the  tunnel  test  section.  The  defined 
in  this  way  is  not  numerically  comparable  to  that  of  Reference  6,  in 
which  experimental  results  were  reported  for  the  same  airfoil  in  a 
nonuniformly  sheared  flow  with  free  boundaries  similar  to  that  of 
Figure  1(c).  The  numerical  difference  arises  because  of  the  difficulty 
of  defining  an  average  U  which  is  equivalent  for  the  two  sheared  flows. 

In  spite  of  the  numerical  differences,  the  shapes  of  the  5/  versus  Qc, 
curves  of  Figure  18  are  similar  to  those  presented  in  Reference  6.  For 
both  sheared  flows,  at  a  -  0,  there  is  a  positive  lift  for  h  >  0  and  a 
negative  lift  for  h  <  0.  The  behavior  of  cj  nc  maximum  lift  appears 
to  be  dependent  on  the  magnitude  of  the  shear.  At  positions  well  below 
the  centerline  of  the  velocity  distribution,  the  stall  is  abrupt.  At 
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positions  near  the  centerline  and  above  the  centerline,  the  onset  of  stall 
becomes  much  more  gradual. 


The  maximum  usable  lift  as  defined  in  Reference  6  is  the  value 

of  cj  at  which  the  derivative  of  the  moment  coefficient  about  the  midchord 

with  respect  to  (  )  ,  i.  e.  ,  —  ,  becomes  zero.  An  apparent 

*  d(cz)R 

correlation  of  this  maximum  usable  lift  with  the  product  of  shear  pa- 
rameter  and  gradient  of  shear  parameter,  (K—),  was  noted  in  Refer- 

dy 

ence  6.  The  present  data  (Figure  18)  showed  no  evidence  of  such  a 
correlation. 

Figure  19  shows  (cj  )*  versus  d  for  all  positions  of  the  model 
that  were  tested  in  the  sheared  flow.  These  data  are  directly  comparable 
to  the  data  of  Figure  6  of  Reference  6  (which  has  been  reproduced  here 
as  Figure  20),  if  the  value  of  r  in  Figure  20  is  taken  to  be  12  inches.  The 
test  configuration  in  the  present  case  differs  from  that  of  Reference  6  in 
the  lower  values  of  shear  and  the  constrained  edges  of  the  nonuniformly 
sheared  regions  of  the  flow.  The  trends  evident  in  Figures  19  and  20  are 
similar.  Note  that  even  in  the  low-shear  tests  (Figure  19),  lift  coefficients, 
(c/ )* ,  as  high  as  2.0  are  obtained.  However,  the  pronounced  effect  of 
different  airfoil  position  near  the  centerline  of  the  nonuniformly  sheared 
free  jet  on  the  shape  of  the  lift  coefficient  curve  (Figure  20,  h/r  -  +  1/8 
and  -1/8)  does  not  occur  with  the  bounded  nonuniform  flow  (Figure  19). 
Moreover,  the  maximum  usable  lift  for  the  bounded  nonuniform  flow  does 
not  show  any  marked  variation  with  airfoil  position  near  the  centerline  of 
the  flow  (Figure  18  does  show  large  variations  in  maximum  usable  lift, 
but  these  variations  correspond  to  relatively  large  changes  in  airfoil 
position.)  This,  again,  is  in  contrast  to  Reference  6,  whf  e  the  maximum 
usable  lift  varied  from  sz  O.r  at  —  =  +  1/8  to  c2^/.f  at  -  -1/8. 

These  large  variations  in  aerodynamic  behavior  with  the  airfoil  at  different 
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positions  near  the  centerline  of  the  free  nonuniform  jet  may  be  a  conse¬ 
quence  of  the  higher  shear  or  of  the  fact  that  free  boundaries  are  present, 
or  a  combination  of  both. 

Figure  21  shows  the  drag  coefficient  (c^  )g  plotted  versus  (  C  t  )*. 
As  mentioned  previously,  (  C ^  )*  includes  the  pressure  drag  only  and  is 
probably  not  very  accurate.  Of  interest  is  the  fact  that  (  )g  is  positive 
in  all  cases.  Sheared  flow  aerodynamic  data  obtained  previously  with 
the  strain-gauge  balance  (Reference  8)  had  resulted  in  negative  drag 
coefficients  for  some  cases  with  small  values  of  (  C/  )g  . 

The  moment  coefficient  about  the  quarterchord,  (cm  /()R  ,  is 
shown  plotted  versus  (  )g  in  Figure  22.  These  data  show  a  relatively 
large  variation  with  differing  wing  positions  and  appear  to  fall  roughly 
into  two  groups,  one  for  data  obtained  with  the  wing  below  the  center- 
line  of  the  sheared  flow  and  one  for  the  data  obtained  with  the  wing  above 
the  centerline.  The  data  obtained  with  h  =  -1/2  inch,  which  was 
closest  to  the  centerline  of  the  experimental  velocity  distribution,  fall 
approximately  midway  between  these  groups.  The  differences  in  the 
groups  of  curves  indicate  variations  in  pressure  distribution  at 
constant  values  of  (  Cj  )R .  Notable  in  Figure  22  is  the  evidence  that  the 
center  of  pressure  moves  continuously  rearward  with  increasing  (  )R 
for  negative  values  of  h  »  while  for  positive  h  ,  it  remains  approx¬ 
imately  constant  or  even  moves  forward  as  (  )g  increases  up  to 
values  near  the  maximum.  Just  below  the  maximum  values  of  (  )R, 

the  center  of  pressure  of  the  ( +h )  group  begins  to  move  rearward 
abruptly  in  a  fashion  similar  to  the  (-/>)  group. 

It  is  interesting  to  compare  pressure  distributions  in  the  regions 
where  the  above  two  groups  show  dissimilar  behavior  in  the  movement 
of  center  of  pressure.  This  has  been  done  in  Figure  23  for  h  -  +3 
and  -3  inches  in  the  sheared  flow  and  also  for  uniform  flow.  The 
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data  shown  were  chosen  to  be  at  values  of  (  ).  as  close  to  each  other 

as  possible  and  in  a  range  of  (  )K  where  separation  of  the  boundary 
layer  should  not  be  too  influential  in  determining  the  pressure  distri¬ 
bution.  Figure  23  shows  that  the  adverse  pressure  gradient  aft  of  the 
suction  peak  on  the  upper  surface  is  increased  by  negative  shear 
(  h  =  -3  inches)  and  decreased  by  positive  shear  ( h  =  +3  inches)  in 
comparison  with  that  which  is  obtained  in  uniform  flow.  This  behavior 
is  consistent  with  the  differences  in  the  shapes  of  the  (C^)^  versus  <x 
curve  at  high  angles  of  attack  for  positive  and  negative  shear.  The 
steeper  pressure  gradient  will  tend  to  promote  separation;  the  lower 
pressure  gradient,  to  delay  separation. 

The  upper  surface  pressure  distributions  of  Figure  23,  which 
are  typical  of  all  values  of  (Cj  )g  up  to  the  maximum,  do  not  explain, 
however,  the  seemingly  haphazard  variation  in  the  maximum  values  of 
(  Ci  )#  illustrated  in  Figure  19.  From  these  upper  surface  pressure 
distributions  it  could  be  concluded  that,  since  separation  is  promoted 
by  negative  shear  and  delayed  by  positive  shear,  the  maximum  values 
of  (  cj  )  should  be  observed  when  the  shear  is  positive  and  a 
progressive  deterioration  in  maximum  ( c/ )-  should  occur  as  the  shear 
becomes  negative.  This  is  certainly  not  the  case  in  Figure  19  where 
the  data  for  h  -  +  3  and  -3  inches  gave  approximately  the  same 
maximum  value  of  (  Cj  )^  ,  although  it  must  be  admitted  (as  mentioned 
in  the  previous  paragraph)  that  the  abrupt  stall  for  negative  shear  and  the 
gradual  stall  for  positive  shear  are  consistent  with  the  upper  surface 
pressure  distributions.  A  partial  explanation  for  this  apparent 
inconsistency  is  that  the  lift  also  depends  on  the  lower  surface  pressures, 
and  as  can  be  seen  in  Figure  23,  the  lift  due  to  the  lower  surface  for 
h  =  -3  is  appreciably  greater  than  that  for  h  =  +3. 

Some  degree  of  correlation  can  be  brought  to  the  data  by  using 
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the  apparent  free-strea  n  dynamic  pressure  of  the  stagnation  streamline 
which  intersects  the  airfoil  as  a  reference  dynamic  pressure.  The  method 
for  obtaining  this  dynamic  pressure  was  discussed  previously.  It  amounts 
to  choosing  the  value  of  dynamic  pressure  which  will  make  the  positive 
peak  of  ( Cp  )*  equal  to  unity.  For  example,  in  Figure  23,  for  h  -  -3 
inches,  all  values  of  (  C?  )4  are  converted  to  {  C p  )5  simply  by  dividing 
(  Cp  by  1. 25.  The  same  factor  is  then  used  to  obtain  (  )s  from  ( )  . 

If  in  addition  to  using  (c^ )  ,  the  angle  of  attack  for  zero  lift  is  subtracted 
from  the  usual  angle  of  attack,  the  resulting  curves  for  the  various  values 
of  h  become  directly  comparable.  The  results  are  plotted  in  Figure 
24(a)  as  (  c*  )5  versus  oi  I  <X  -  &Ci  a-  A  similar  presentation  of  ( )R 
is  shown  in  Figure  24(b)  for  comparison.  Note  that  individual  sets  of 
(  Cj,  )s  data  for  a  given  airfoil  midchord  position  show  a  greater  degree 
of  irregularity  (they  do  not  form  as  smooth  a  curve)  than  the 
corresponding  (  cx  )/r  data.  This  can  be  attributed  to  errors  in  determining 

the  free -stream  dynamic  pressure  of  the  stagnation  streamline  from  the 
experimental  pressure  distributions.  Unfortunately,  there  was  a 
greater  concentration  of  static-pressure  taps  on  the  upper  leading- 
edge  surface  of  the  airfoil  than  on  the  lower  leading  edge.  Consequently, 
in  many  cases,  there  was  difficulty  in  determining  accurately  the 
maximum  positive  value  of  (  Qp  )^ » and  such  inaccuracies  are  reflected 
directly  in  (  c»  )  . 

Inspection  of  the  (  )  data  (Figure  24(b))  provides  several 

features  which  are  worthy  of  note.  First,  all  of  the  experimental  data 
obtained  in  the  sheared  flow  for  values  of  ac  greater  than  zero  lie  above 
the  control  data  for  the  same  airfoil  in  uniform  flow  at  approximately 
the  same  Reynolds  number.  Second,  the  slopes  (  )  of  a  number 

of  the  sets  of  sheared-flow  data  increase  as  angle  of  attack  increases, 
becoming  greater  than  that  predicted  by  inviscid  uniform  flow  theory. 

This  latter  observation  cannot  be  accounted  for  on  the  basis  of  viscous 
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effects  alone,  since  such  effects  on  an  airfoil  always  tend  to  result  in 
reduced  lift  compared  to  that  predicted  by  inviscid  theory. 

Inspection  of  Figure  24(a)  provides  a  somewhat  different,  and 
more  enlightening, picture.  In  this  case,  all  of  the  experimental  data 
at  positive  values  of  (X  lie  on  or  below  the  inviscid  uniform  flow 
theoretical  curve  and  appear  to  fall  into  two  approximate  groupings,  one 
for  positions  below  the  sheared  flow  centerline  and  one  for  positions 
near  and  above  the  sheared  flow  centerline.  These  groupings  are 
similar  to  those  noted  previously  in  the  (  Cm  versus  ( c/  data 
(Figure  22).  Such  correspondence  is  not  particularly  surprising,  since 
plotting  versus  cz  effectively  cancels  any  error  which  arises  through 
the  use  of  incorrect  or  inconsistent  values  of  reference  dynamic  pressure 
and  completely  eliminates  the  angle  of  attack  from  the  presentation. 

The  differences  between  the  two  groups  of  data  in  Figure  24(a)  will  be 
shown  to  be  consistent  with  differences  in  the  observed  pressure 
•  jtributions. 

Figure  25  shows  the  pressure  distribution  (  CF  )9  on  the  airfoil 
for  different  wing  positions.  The  distributions  shown  are  for  nearly 
equal  values  of  a,  with  cc  small  enough  that  boundary-layer  separation 
should  not  affect  the  results  to  any  marked  degree.  Note  that  the 
pressure  gradient  after  the  suction  peak  is  less  steep  for  positive 
shear  (  h  =  +3  inches)  than  for  negative  shear  {h  -  -3  inches)  and 
also  that  the  gradient  shown  by  the  uniform  flow  control  data  is 
approximately  the  same  as  that  for  the  negative  shear  data.  The 
magnitude  of  (  C p  )3  at  the  suction  peak  is  influenced  by  the  existence 
of  a  leading-edge  bubble.  This  bubble  occurred  both  with  uniform 
flow  and  with  h  =  -3  inches, even  at  low  angles  of  attack.  At  h  -  +3 
inches,  however,  the  bubble  did  not  appear  until  the  angle  of  attack 
approached  10  degrees.  Hence,  comparison  of  the  magnitude  of  the 
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suction  peaks  is  not  of  any  particular  value.  Comparison  of  the  pressure 
gradients  after  the  suction  peak  is  useful  because  the  uniform  flow 
control  experiments  indicated  that  these  gradients  are  not  particularly 
sensitive  to  the  existence  of  the  bubble.  (See* for  example,  Figures  10, 

1 1 ,  and  12 . ) 

Similar  pressure  distributions  are  shown  in  Figure  26  for  a 
higher  value  of  a.  The  same  trends  are  evident  here  as  were  displayed 
in  Figure  25.  Hence,  it  appears  that  the  shapes  of  the  (  )  curves 

near  stall  are  consistent  with  the  pressure  distributions.  These 
distributions  indicate  that  separation  should  be  delayed  by  positive 
shear  and  promoted  by  negative  shear.  For  the  airfoil  used  in  these 
tests,  the  sharp  suction  peaks  and  large  adverse  pressure  gradients 
generated  in  the  uniform  flow  and  the  negatively  sheared  flow  cause 
the  formation  of  a  leading-edge  bubble  which  reduces  the  lift  to  values 
less  than  those  predicted  by  uniform  flow  theory  even  at  relatively  low 
angles  of  attack.  The  favorable  effect  of  positive  shear  dr-)\y  s  the 
formation  of  the  leading-edge  bubble  and  also  delays  th^  onset  of 
trailing-edge  separation, resulting  in  a  closer  correlation  with  uniform- 
flow,  inviscid  theory  in  Figure  24(a). 

From  the  above  discussion,  it  would  appear  that  a  large  portion 
of  the  confusion  regarding  the  effects  of  nonuniform  shear  on  the  lift 
of  an  airfoil  can  be  eliminated  if  a  proper  reference  velocity  is  used  to 
reduce  the  data  to  nondimensional  form.  This  velocity  is  the  free- 
stream  velocity  which  corresponds  to  the  stagnation  streamline  of  the 
airfoil. 


Incidental  to  the  determination  of  the  free -stream  dynamic 
pressure  of  the  stagnation  streamline,  displacement  of  the  stagnation 
streamlines  at  the  various  angles  of  attack  for  each  value  of  K  tested 
was  also  determined.  (See  sketch  on  following  page. ) 
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Stagnation  streamline  displacement,  &YS  ,  is  the  vertical  distance  above 
the  airfoil  midchord  reference  height,  h  ,  to  the  location  in  the  free-stream 
velocity  profile  where  U(y)  £  Us  .  The  ability  to  trace  the  approximate 
location  from  which  the  stagnation  streamline  originated  through 
simple  pressure  measurements  on  a  body  (ignoring  viscous  losses  due 
to  mixing)  is  a  unique  and  useful  feature  of  two- dimensional  sheared 
flows. 


The  resultant  data  for  h  -  -3,  -2,  -1,  +1,  and  +3  inches  are 
plotted  in  Figures  27  through  31.  Again,  these  data  should  be  viewed 
with  some  caution,  as  they  were  derived  from  the  stagnation  point 
static  pressure  and  the  errors  may  be  large.  Bars  are  drawn  through 
some  of  the  particularly  questionable  points  to  indicate  the  possible 
degree  of  error  arising  from  estimating  the  maximum  positive  value  of 
(  Cp  )p  .  In  addition  to  errors  due  directly  to  the  difficulty  in  obtaining 
the  correct  value  of  [Cp  )p  at  the  stagnation  point,  there  was  further 
magnification  of  these  errors  when  the  airfoil  was  near  the  centerline 
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of  the  sheared  flow.  The  rounded  velocity  profile  in  this  region  made 
small  errors  in  (Cp  )*  at  the  stagnation  point  correspond  to  large  errors 
in  AYS  .  The  values  of  AYS  derived  experimentally  for  h  =  +1/2,  0, 

-1/2  were  particularly  suspect  because  of  this  and,  hence,  are  not 
presented.  The  data  for  the  wing  well  below  the  centerline  of  the  sheared 
flow  behave  in  a  manner  which  appears  to  be  consistent  with  an  increase 
in  circulation  around  the  airfoil  as  angle  of  attack  is  increased.  The  data 
obtained  with  the  airfoil  above  the  centerline  vary  in  a  manner  different 
from  the  data  obtained  below  the  centerline.  No  explanation  for  the 
differences  in  behavior  of  the  apparent  deflection  of  the  stagnation  flow 
streamline  is  advanced  at  present.  Further  computations  using  the  non- 
uniform  sheared  flow  theory  may  help  to  explain  these  differences  in 
stagnation  streamline  deflection. 

From  the  experimental  data  presented  thus  far,  several 
important  conclusions  can  be  stated  with  respect  to  two-dimensional 
aerodynamic  characteristics  in  the  specific  nonuniformly  sheared  flow 
considered  during  the  present  research; 

1.  The  addition  of  turbulence  to  the  flow  through  the  use 
of  a  screen  to  generate  the  sheared  flow  is  not  the 
source  of  the  observed  differences  between  stalling 
characteristics  of  airfoils  in  sheared  and  in  unsheared  flows. 

2.  Much  of  the  apparent  difference  in  aerodynamic  behavior 
below  stall  resulting  from  differing  vertical  positions 

in  the  sheared  flow  can  be  reconciled  if  the  free-stream 
velocity  of  the  stagnation  streamline  which  intersects 
the  airfoil  is  used  as  a  reference  velocity  in  calculating 
lift  coefficients;  see  Figure  24(a).  It  is  worth  noting 
that,  even  in  the  relatively  low  shear  of  the  present 
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work,  the  dynamic  pressure  of  the  stagnation  streamline 
differed  from  the  dynamic  pressure  referenced  to  the 
airfoil  midchord  by  more  than  50  percent  in  some  instances. 

3.  There  are  differences  in  the  approach  to  complete  stall 
of  the  airfoil  depending  upon  the  vertical  location  of  the 
airfoil  in  the  sheared  flow.  These  differences  appear  to  be 
associated  with  the  change  in  pressure  distribution  on  the 
upper  surface  of  the  airfoil.  Positive  shear  at  the  airfoil 
position  reduces  the  pressure  gradient  on  the  upper  surface 
and  hence  delays  the  occurrence  of  separation,  while  the 
converse  is  true  of  negative  shear. 

These  observations  suggest  that  it  is  probably  not  necessary  to 
consider  the  direct  interaction  between  shear  in  the  free  stream  and  the 
behavior  of  the  boundary  layer  when  seeking  the  stalling  characteristics  of 
a  wing  in  a  sheared  flow.  If  the  pressure  distribution  on  a  wing  in  an 
inviscid  sheared  flow  is  predictable,  then  the  prediction  of  stalling 
characteristics  in  a  real  sheared  flow  should  be  equivalent  to  the 
prediction  of  stalling  characteristics  of  a  wing  in  a  viscous, uniform 
flow,  given  the  inviscid, uniform-flow  pressure  distribution.  This 
hypothesis  will  be  further  substantiated  by  comparison  of  theoretical 
and  experimental  pressure  distributions  in  the  sheared  flow  in  subsequent 
paragraphs. 
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COMPUTED  RESULTS  AND  COMPARISON  WITH  EXPERIMENTAL  RESULTS 


Calculations  using  the  developed  computer  program  for  a  non- 
uniformly  sheared  flow  were  performed  corresponding  to  two  vertical 
locations  of  the  wing  with  respect  to  the  flow  centerline,  h/c  -  -0.33 
and  h/c  -  +0.  50.  The  free-stream  velocity  profile  used  in  the  numerical 
calculations  was  such  that  the  shear  parameter,  K,  varied  from  i  1  u^ar 
the  tunnel  centerline  to±0.  33  adjacent  to  the  wind-tunnel  walls.  At  K/c  - 
-0.  33,  pressure  distributions  were  computed  for  cl  =  7.  1  degrees  and 
oc=  9.  5  degrees.  At  Vy£  =  +0.  5,  the  pressure  distribution  at0L=  8.  2 
degrees  was  computed. 

Figure  32  is  a  comparison  of  computed  and  experimental  values 
of  (  Cp  versus  z/c  for  h  =  -2  inches  (corresponding  to  ^  =  -0.  33), 

<2=  7.  1  degrees;  Figure  33, for  h  =  - 2  inches,  cc  ~  9.  5  degrees 
(experimental  data  for  (X  =  9.  6  degrees);  and  Figure  34,for  h  =  +3 
inches  (  h/c  =  0.  5),  (X  -  8.  2  degrees.  It  can  be  seen  that  in  all  three 
cases  the  comparison  is  as  good  as  in  the  corresponding  uniform-flow 
case  between  the  results  of  inviscid  theory  and  the  experimental  data  in 
a  low  Reynolds  number  turbulent  flow  (see  Figures  10  through  12).  In 
particular,  the  upper-surface  pressure  gradient  is  quite  accurately 
predicted,  although,  as  for  the  uniform-flow  cases,  the  leading  edge 
suction  peak  is  overpredicted,  more  so  as  the  peak  ( Cp)#  becomes 
more  negative. 

In  Figure  34,  for  h  -  +3  inches,  the  uniform-flow  pressure 
distribution  calculated  from  the  computer  program  for  the  airfoil  at 
h/c  =  +0.  5  is  also  plotted  for  comparison.  It  is  clear  that,in  agree¬ 
ment  with  experiment,  the  predicted  pressure  gradient  just  aft  of  the 
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leading-edge  suction  peak  is  less  adverse  for  the  h  -  +3  inches,  nonuni- 
formly  sheared-flow  configuration  than  for  the  uniform-flow  configuration. 

Computed  (Cj  )#  and  (c*,0/y)*  for  these  three  cases  are  plotted 
on  Figures  35  and  36,  respectively,  where  comparisons  are  made  with 
the  appropriate  experimental  results.  Agreement  is  satisfactory  except, 
perhaps,  for  h  -  - 2  inches,  or  -  9.  5  degrees  {(Cj  )*  =  1. 72),  where  experi¬ 
mentally  the  center  of  pressure  is  further  aft  than  predicted.  Computed 
values  of  )f,  including  the  uniform  flow  configuration,  vary  from 
-0.010  to  -0.016.  This  reflects  the  slight  overprediction  of  leading-edge 
peak  suction  for  uniform  flow,  which  undoubtedly  carries  over  into  the 
nonuniformly  sheared-flow  calculations. 

Computed  displacements  of  the  airfoil  stagnation  streamlines  are 
plotted  in  Figures  28  and  31.  Relatively  good  agreement  with  the  experi¬ 
mentally  derived  displacements  is  obtained  for  h  -  -2  inches  (Figure  28). 
Predicted  displacements  in  this  case  are  25  percent  to  35  percent  too  high. 
Poorer  agreement  is  shown  for  the  h  =  +3  inches,  Cc  -  8.  2  degrees  case 
(Figure  31).  If  there  were  perfect  symmetry  in  the  experimental  nonuniform 
velocity  profile,  negative  stagnation  streamline  deflections  for  the  airfoil 
at  h  -  -3  inches  at  positive  angles  of  attack  would  correspond  to  positive 
streamline  deflections  with  the  airfoil  at  h  -  +3  inches  at  negative  angles 
of  attack.  The  experimental  stagnation  streamline  displacement  data  for 
h  -  -3  inches,  with  signs  of  as  and  &  Ys  reversed,  is  also  plotted  in  Figure  31. 
The  comparison  of  computed  and  experimental  stagnation  streamline  de¬ 
flection  data  is  somewhat  improved  in  this  case.  (Note  that  the  corre¬ 
sponding  pressure  distribution  for  h  =  -3  inches,  ot>  =  -8.2  degrees  is  plotted 
in  Figure  34  for  comparison  with  h  -  +3  inches,  a  -  8.  2  degrees  data. )  In 
view  of  the  uncertainties  involved  in  the  determination  of  experimental 
stagnation  streamline  deflections,  the  comparison  of  computed  and 
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experimental  results  in  this  case  is  quite  good. 


DISCUSSION  OF  COMPUTED  RESULTS 

It  was  concluded  on  the  basis  of  the  experimental  data  that  the 
variation  in  stall  behavior  with  vertical  location  of  the  airfoil  in  the 
nonuniformly  sheared  flow  was  associated  with  the  corresponding 
variation  in  the  adverse  pressure  gradient  aft  of  the  suction  peak  on  the 
airfoil  upper  surface.  For  the  airfoil  above  the  centerline  of  the  non¬ 
uniformly  sheared  flow,  the  onset  of  trailing -edge  stall  was  delayed 
somewhat,  and  the  occurrence  of  complete  separation  was  delayed  to 
somewhat  higher  angles  of  attack  than  in  the  case  of  a  directly  comparable 
uniform  flow  (at  the  same  Reynolds  number  and  with  similar  turbulence 
level).  Related  to  this  behavior  was  a  less  adverse  pressure  gradient 
on  the  upper  surface  of  the  airfoil  in  the  nonuniformly  sheared  flow  than 
in  the  uniform  flow.  Likewise,  for  the  airfoil  below  the  centerline  of 
the  nonuniformly  sheared  flow,  stall  was  correspondingly  promoted,  and 
in  this  case,  the  upper  surface  pressure  gradient  was  more  adverse  than 
for  the  comparable  uniform  flow. 

The  question  remains:  are  the  experimentally  demonstrated 
variations  of  adverse  pressure  gradient  with  change  in  vertical  height 
associated  with  an  interaction  of  the  sheared  flow  and  the  boundary 
layer,  or  are  they  attributable  solely  to  the  modification  of  the 
pressure  gradient  by  the  sheared  flow?  (Note  that  the  laminar 
separation  bubble  is  an  example  of  an  interaction  of  the  boundary  layer 
and  the  outer  flow.  )  The  good  agreement  of  the  computed  and  experi¬ 
mental  pressure  distributions  indicates  that,  in  fact,  the  altered  pressure 
distributions  are  primarily  a  result  of  inviscid  action  of  the  sheared 
flow. 


The  flow  model  upon  which  the  digital  computer  program  is 
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based  involves  several  idealizations  which  could  affect  the  computed 
results  to  some  degree.  The  two  most  important  are  probably: 

1.  The  flow  is  assumed  to  be  inviscid. 

2.  The  shear  screen  is  assumed  to  be  located  a  great 
distance  upstream  relative  to  the  model  chord  length. 

3.  Flow  disturbances  are  assumed  localized  in  the 
immediate  vicinity  of  the  airfoil. 

The  assumption  of  an  inviscid  fluid  means  that  both  viscous 
mixing  and  the  displacement  effect  of  the  boundary  layer  are  neglected. 

If  the  disturbance  arising  from  the  airfoil  is  localized  in  the  flow  within 
two  or  three  chord  lengths  of  the  airfoil,  the  effects  of  neglecting  viscous 
mixing  in  the  theoretical  model  are  probably  small.  Indications  from 
the  limited  calculations  done  thus  far  are  that  the  airfoil  disturbance  is 
felt  farther  upstream  and  downstream  for  nonuniformly  sheared  flow  than 
for  uniform  flow,  particularly  for  the  h  =  +3  inches,  <x  =  8.  2  degrees 
configuration. 

The  effect  of  assuming  that  the  shear  screen  is  very  far  upstream 
of  the  model  is  difficult  to  ascertain.  In  actuality,  the  shear  screen  was 
located  about  reven  airfoil  chord  lengths  upstream  of  the  airfoil  midchord. 

If  the  airfoil  flow  disturbance  did  extend  this  far  upstream,  then  there 
could  be  an  interaction  of  some  sort.  Theoretically,  it  would  have  been 
possible  to  include  the  shear  screen  in  the  flow  model;  however,  limitations 
of  computer  capacity  precluded  this  for  the  present. 

The  assumption  that  flow  disturbances  associated  with  the  airfoil 
do  not  propagate  indefinitely  far  upstream  or  downstream  of  the  airfoil  in 
the  wind  tunnel  is  incorporated  in  the  computer  program  by  fixing  g  (*)  — 
the  displacement  of  the  free  streamline  at  the  apex  of  the  velocity  distri¬ 
bution  —  at  its  free-stream  value  (zero  in  the  present  nonuniformly 


39 


sheared  flow)  for  z  <-.£5  and  z>  >  Ds  ,  and  by  fixing  the  vortex  strength 
distribution  representing  the  wind-tunnel  wall  at  the  free -stream  value 
for  z  <-^and  z>  D# .  For  the  calculations  at  h  ~  -2  inches,  the  values 
(ps)A=  (o„)/c  =4  were  chosen.  For  the  h  =  +3  inches  configuration, 
calculations  were  performed  for  {Ps)/c  =  ( Dw)jc  -  6  and  for  (p 't)Jo  = 

{P*)/c  -  8.  In  principle,  the  extent  of  the  flow  disturbance  upstream 
and  downstream  of  the  airfoil  could  be  computed  by  arbitrarily  choosing 
values  of  Pw  and  Ps  ,  obtaining  a  solution,  and  determining  whether,  in 
fact,  the  vortex  strengths  for  the  wind-tunnel  wall  at  z  -  *  Pw  and  gfe) 
at  X  =  *  Ps  have  attained  their  free -stream  values.  If  not,  Ps  and  D# 
could  be  made  larger  and  the  calculations  repeated  until  values  of  Ds 
and  Dh  were  found  for  which  ^  (*  As  )»  (  *  D#  ),  and  yz  (  *  0^)  had 

attained  their  free-stream  values.  This  was  not  practicable  because  of 
the  limitations  of  computer  storage  capacity  and  also  the  practical 
limitation  on  the  maximum  size  of  matrix  that  can  be  inverted 
numerically  with  acceptable  accuracy.  It  was  found  for  the  h  =  +3  inches 
calculations  that  as  ^  and  P#  were  increased  from  six  to  eight  chord 
lengths,  the  values  of  g{±Ds),  (±0^),  and  showed  even  less 

agreement  with  free-stream  values.  This  behavior,  at  least  ins*-'  * 

Yf  (  and  yz  (  *  P^)  are  concerned,  was  also  noted  in  sevei  e 

uniform  flow  cases  that  were  computed  with  the  program.  For  this 
reason,  it  is  believed  that  this  result  may  partly  be  associated  with  the 
numerical  techniques  used  in  the  computer  program.  In  particular,  the 
small  numerical  inaccuracies  associated  with  the  smoothing  procedure 
might  have  been  a  majcr  contributing  factor.  In  any  case,  these  noted 
discrepancies  were  sma11  in  the  present  calculations. 

The  comparisons  of  theoretical  and  experimental  data, 
particularly  the  pressure  distributions,  indicate  that  the  theory 
successfully  predicts  the  important  features  of  the  (unseparated  flow) 
airfoil  aerodynamic  characteristics  in  the  nonuniformly  sheared  flow 
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treated  during  the  present  research.  It  appears  that  the  stall  behavior 
can  be  accounted  for  on  the  basis  of  the  inviscid  flow  pressure  distri¬ 
bution,  and  that  there  are  no  unusually  large  viscous  interaction  effects 
between  the  sheared  outer  flow  and  the  boundary  layer.  The  analytical 
technique,  therefore,  should  prove  to  be  an  extremely  useful  research 
tool  for  further  investigation  of  the  stall  behavior  of  airfoils  in  nonuniformly 
sheared  flows. 

It  would  be  most  desirable  to  extend  the  computed  results  for  the 
present  velocity  profile  over  a  greater  range  of  angles  of  attack  and  for 
additional  />/c  values.  This  would  permit  a  theoretical  check  of  the  rela- 

fa 

tively  good  correlation  of  the  experimental  data  for  -0.  5  £  —  <  0.  5  in 
terms  of  (c^)s  versus  a  .  In  particular,  it  would  be  of  great  interest  to 
determine  how  accurately  the  theory  predicts  ( Ac >)-  at  oc  -  0.  The  cal¬ 
culations  should  also  be  extended  to  higher  values  of  the  shear  parameter 
K  than  those  considered  in  the  current  work. 

The  experimental  airfoil  data  of  Reference  6  were  obtained  in 
a  two-dimensional,  nonuniformly  sheared  flow  with  nearly  four  times 
the  shear  near  the  flow  centerline  as  for  the  flow  of  the  present  tests. 

The  sheared  flow  computer  program  should  be  modified  and  extended 
to  the  nonuniform  velocity  profile  of  the  Reference  6  data,  and  pressure 
distributions  should  be  computed  for  a  number  of  the  configurations  of 
the  tests  of  Reference  6.  The  shear  in  the  Reference  6  velocity  profile 
is  probably  as  large  or  larger  than  the  maximum  obtainable  in  a 
propeller  slipstream.  If,  as  expected,  tne  theory  retains  its  accuracy 
for  this  level  of  shear,  its  range  of  validity  should  be  adequate  to  cover 
all  foreseeable  practical  cases. 


PROSPECTS  FOR  ACCOUNTING  FOR  THREE-DIMENSIONAL  EFFECTS 


An  important  application  of  the  results  obtained  from  this 
continuing  research  program  is,  of  course,  to  the  aerodynamics  of 
wings  immersed  in  propeller  slipstreams.  The  engineering  prediction 
of  the  lift  generated  by  a  wing  interacting  with  a  propeller  slipstream 
at  low  forward  speeds  has  been  a  difficult  problem.  This  problem  is 
of  particular  importance  for  VTOL  and  STOL  aircraft  in  which  either 
the  propeller  thrust  or  the  lift  generated  by  deflection  of  propeller 
slipstreams  is  a  substantial  portion  of  the  take-off  lift  of  the  aircraft. 

There  are  at  least  two  possible  approaches  to  the  analytical 
prediction  of  lift  on  a  finite  wing  in  a  propeller  slipstream  : 

1.  A  direct  approach  based  on  a  vortex  sheet  flow  model 
similar  to  that  of  the  analytical  technique  used  in  the 
present  study. 

2.  A  lifting-line  approach  based  on  the  theory  for  three- 
dimensional,  nonuniformly  sheared  flows  proposed  by 
von  Karman  and  Tsien  (Reference  13). 

Clearly,  the  first  approach  would  result  in  a  computer  program  of 
formidable  complexity  even  in  comparison  with  the  two-dimensional 
flow  computer  program  developed  during  the  present  research.  The 
more  attractive  approach  from  a  practical  computational  point  of  view 
is  the  lifting -line  approximation.  However,  even  the  von  Karman- 
Tsien  theory  cannot  be  applied  readily  to  the  calculation  of  finite-wing 
aerodynamic  characteristics  because  the  theoretical  expression  for 
the  induced  angle  of  attack  at  each  wing  section  must  be  rederived  for 
each  nonuniformly  sheared  flow  velocity  distribution.  The  derivation 
involves  the  solution  of  two  differential  equations  and  an  integral 
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equation;  finding  these  solutions  undoubtedly  will  not  be  a  simple  process. 
In  addition,  the  von  Karman-Tsien  theory  is  a  small  disturbance  theory, 
and  its  extension  to  high  wing  loadings  near  maximum  lift  can  be  justified 
only  on  the  basis  of  a  comparison  of  the  results  of  analysis  and  experiment. 

In  any  case,  two-dimensional  sectional  lift  data  derived  from 
two-dimensional, nonuniformly  sheared  flows  with  the  same  velocity 
profiles  as  those  in  sections  vertical  to  the  plane  of  the  finite  wing  are 
needed  in  order  to  apply  the  lifting-line  approximation.  If  experimental 
data  were  to  be  obtained  for  this  purpose,  a  difficult  and  expensive  test 
program  would  be  required.  An  accurate  two-dimensional  airfoil  theory 
for  nonuniformly  sheared  flow  would  provide  the  same  sort  of  data,  and 
at  considerably  lees  expense. 

It  is  believed  that  effort  should  be  devoted  to  exploitation  of  the 
von  Karman-Tsien  nonuniform  flow  lifting-line  theory.  If  this  effort 
were  successful,  then  the  two-  and  three-dimensional  experimental 
flow  data  developed  and  presented  in  References  6  and  8  would  permit 
an  evaluation  of  the  usefulness  of  such  a  modified  lifting-line  theory  as 
an  engineering  tool.  Additional  two-dimensional  sectional  data  for  the 
nonuniform  flow  corresponding  to  the  vertical  velocity  profiles  of  the 
three-dimensional  nonuniform  flow  treated  in  Reference  8  will 
undoubtedly  be  required.  These  data  could  be  generated  from  the  two- 
dimensional  theory  of  this  report. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


On  the  basis  of  the  work  done  and  the  results  obtained,  the 
following  conclusions  can  be  drawn  relative  to  the  aerodynamic  prop¬ 
erties  of  a  Joukowski  17- percent-thick  symmetric  airfoil  in  the  two- 
dimensional  nonuniformly  sheared  flow  considered  in  this  program 
(uniform  positive  shear  above  the  flow  centerline ;  uniform  negative 
shear  below  the  flow  centerline). 

1.  The  addition  of  turbulence  to  the  flow  through  the  use  of  a 

screen  to  generate  the  sheared  flow  is  not  the  source  of  the  observed 
differences  between  stalling  characteristics  of  airfoils  in  sheared  and 
in  unsheared  flows. 

2.  The  airfoil's  stall  behavior  varies  with  vertical  position  of 
the  airfoil  in  the  flow.  For  airfoil  positions  above  the  flow  centerline, 
(positive  shear)  stall  is  delayed;  for  the  airfoil  below  the  flow  center  - 
line  (negative  shear),  stall  is  promoted. 

3 .  The  stall  behavior  is  consistent  with  the  experimentally 
observed  variations  in  the  adverse  pressure  gradients  on  the  upper 
surface  of  the  airfoil.  Smaller  adverse  pressure  gradients  were  observed 
for  the  airfoil  positioned  above  the  flow  centerline  than  for  the  airfoil 
below  the  flow  centerline. 

4  .  Calculated  pressure  distributions  based  on  an  inviscid  theory 
for  an  airfoil  in  a  nonuniformly  sheared  flow  were  in  good  agreement 
with  experimental  pressure  distributions,  particularly  with  respect  to 
the  upper  surface  pressure  gradients.  It  appears,  therefore,  that  the 
experimentally  observed  behavior  of  airfoils  near  maximum  lift  in  a 


44 


nonuniformly  sheared  flow  can  be  attributed  primarily  to  the  effects  of 
the  inviscid  pressure  distribution  on  the  upper  surface  boundary  layer, 
as  is  the  case  for  airfoils  in  uniform  flows. 

5.  A  marked  degree  of  correlation  of  the  data  is  obtained  for  differ¬ 
ing  vertical  positions  in  the  flow  if  section  lift  coefficients  are  referenced 
to  the  airfoil's  stagnation  streamline  dynamic  pressure,  and  angle  of 
attack  is  measured  from  the  zero -lift  angle  of  attack.  The  two  primary 
factors  in  determining  the  observed  two-dimensional  aerodynamic 
behavior  of  airfoils  in  nonuniformly  sheared  flows,  then,  appear  to  be  the 
variation  of  stagnation  streamline  dynamic  pressure  with  angle  of  attack 
for  a  given  wing  position  in  the  sheared  flow  and  the  inviscid  influence 
of  the  nonuniformly  sheared  flow  in  determining  the  pressure  distribution 
on  the  airfoil  upper  surface.  The  former  factor  determines  the  apparent 
lift-curve  slope  well  below  stall, and  the  latter  determines  the  stalling 
behavior.  Both  factors  appear  to  depend  on  the  overall  distribution  of 
shear  in  the  flow, 

6  .  For  the  particular  nonuniformly  sheared  flow  considered  in 
this  report,  there  was  no  evidence  of  unusual  variations  in  maximum 
lift  obtained  as  airfoil  vertical  position  was  varied  in  the  vicinity  of  the 
flow  centerline.  (This  is  in  contrast  to  the  results  reported  in  Reference 
6  for  a  nonuniformly  sheared  flow  with  larger  shear  and  unconstrained 
flow  boundaries.  ) 

On  the  basis  of  the  results  obtained,  the  following  recommendations 
are  made. 

1 .  Supplemental  airfoil  pressure  distribution  data  should  be 
obtained  experimentally  for  the  particular  two-dimensional,  nonuniformly 
sheared  flow  in  which  the  force  and  moment  data  of  Reference  6  were 
obtained.  Of  special  interest  is  whether  correlation  of  (  CL  )s  versus  Sc 
with  h  similar  to  that  for  the  data  presented  in  this  report  is  obtained. 
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2,  The  computer  program  as  it  is  now  constituted  (constrained 
nonuniformly  sheared  flow  of  Figure  1(b))  should  be  used  for  calculations 
of  airfoil  pressure  distributions  for  additional  values  of  angle  of  attack 
and  vertical  position  of  the  airfoil  in  the  flow.  The  calculations  should 
also  be  extended  to  higher  values  of  the  shear  parameter  K  than  those 
considered  thus  far. 

3.  The  theory  and  computer  program  should  be  extended  to  the 
particular  two  -  dimensional,  nonuniformly  sheared  flow  for  which  the 
data  of  Reference  6  were  obtained.  Calculations  of  airfoil  pressure 
distributions  should  then  be  made  for  a  range  of  airfoil  vertical  positions 
in  the  flow,  and  for  a  range  of  angles  of  attack. 

It  should  be  possible,  from  the  results  of  these  calculations,  and 
from  those  of  (2)  above,  to  determine  the  mechanism  responsible  for  the 
maximum  lift  behavior  reported  in  Reference  6. 
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PRESSURE  COEFFICIENT 


AIRFOIL  COORDINATE,  %/c  ,  FRACTION  OF  CHORD 


Figure  3.  PRESSURE  COEFFICIENT,  CP  ,  VERSUS  FRACTION  OF  CHORD, X/c, 
FOR  UNIFORM  FLOW  FROM  POTENTIAL  FLOW  THEORY,  AND  FROM 
COMPUTER  PROGRAM  WITHOUT  SMOOTH  I  NGj  a  =  10.0  DEGREES. 
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PRESSURE  COEFFICIENT 


0  0.2  0.4  0.6  0.6  1.0 

AIRFOIL  COORDINATE,*/;,  FRACTION  OF  CHORD 


Figure  4.  PRESSURE  COEFFICIENT,  Cp  ,  VERSUS  FRACTION  OF  CHORD,  x/c, 

FOR  UNIFORM  FLOW  FROM  POTENTIAL  FLOW  THEORY,  ANQ  FROM 
COMPUTER  PROGRAM  WITH  SMOOTHING  (*  =  0.I0);<*  =  10.0  DEGREES. 
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Figure  5.  VIEW  OF  AIRFOIL  MODEL  INSTALLED  IN  WIND-TUNNEL  TEST  SECTION. 
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Figure  6.  SHEAR  SCREEN  ASSEMBLY. 


SYMBOL  SPAMWISE  LOCATION 

O  TUNNEL  CENTERLINE 

A  0.75  inch  LEFT  OF  £  (LOOKING  DOWNSTREAM) 


0  20  40  60  80  100  120  140 

VELOCITY, U(y) ,  f.p.t. 

Figure  7.  MEASURED  FREE-STREAM  VELOCITY  DISTRIBUTION  IN  NONUNIFORMLY 
SHEARED  FLOW  AT  MODEL  MIDCHORD  STATION. 
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VERTICAL  HEIGHT  IN  WIND  TUNNEL  FROM  CENTERLINE,  u  ,  inches 


SYMBOL  SPANWISE  LOCATION 

O  TUNNEL  CENTERLINE 

A  0.75  inch  LEFT  OF  {  (LOOKING  DOWNSTREAM) 

□  2.00  inches  LEFT  OF  {  (LOOKING  DOWNSTREAM) 


FREE-STREAM  FLOW  ANGULARITY  REFERENCED 
TO  TEST  SECTION  FLOOR,  degrees 

Figure  8.  MEASURED  FREE-STREAM  FLOW  ANGULARITY  IN  NONUNIFORMLY 
SHEARED  FLOW  AT  MODEL  MIDCHORD  STATION. 


53 


SECTION  U  FT  COEFFICIENT 


ANGLE  OF  ATTACK,  a,  degrees 


Figure  9.  SECTION  LIFT  COEFFICIENT, 4  ,  VERSUS  ANGLE  OF  ATTACK,  a, 
FOR  UNIFORM  FLOW. 
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Figure  10.  PRESSURE  COEFFICIENT,  Cp  ,  VERSUS  FRACTION  OF  CHORD,  */c  , 
UNIFORM  FLOW,  a=  0  DEGREES. 
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Figure  II.  PRESSURE  COEFFICIENT,  Cp  ,  VERSUS  FRACTION  OF  CHORD,  x/c, 
UNIFORM  FLOW,  Or  =  7.5  DEGREES. 
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PRESSURE  COEFFICIENT 


Figure  12. 


PRESSURE  COEFFICIENT,  CF  ,  VERSUS  FRACTION  OF  CHORD,  X/C  , 
UNIFORM  FLOW,  a  =  10.0  DEGREES. 
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PRESSURE  COEFFICIENT 


Figure  13.  PRESSURE  COEFFICIENT,  Cp,  VERSUS  FRACTION  OF  CHORD,  x/c  , 
UNIFORM  FLOW,  LOW  SOLIDITY  TURBULENCE  SCREEN, 

OC=  13.0  DEGREES,  00  =  66  FT./SEC. 
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PRESSURE  COEFFICIENT 
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AIRFOIL  COORDINATE,  %/c  »  FRACTION  OF  CHORD 


Figure  14.  PRESSURE  COEFFICIENT, Cp  ,  VERSUS  FRACTION  OF  CHORD,  */<5, 
UNIFORM  FLOW,  LOW  SOLIDITY  TURBULENCE  SCREEN, 
a  =  16.5  DEGREES,  U0  =  66  FT./SEC. 
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PRESSURE  COEFFICIENT 


Figure  15.  PRESSURE  COEFFICIENT,  C.  ,  VERSUS  FRACTION  OF  CHORD,  */c  , 
UNIFORM  FLOW,  LOW  SOLIDITY  TURBULENCE  SCREEN, 
a=  16.7  DEGREES,  U0  =  66  FT./SEC. 
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SECTION  QUARTERCHORD  MOMENT  COEFFICIENT, 


-6-4-2  0  2  4  6  8  10  12  14  16  18 


ANGLE  OF  ATTACK,  <2,  degrees 


Figure  16.  SECTION  QUARTERCHORD  MOMENT  COEFFICIENT, 

VERSUS  ANGLE  OF  ATTACK,  a,  IN  UNIFORM  FLOW 
WITH  LOW  SOLIDITY  TURBULENCE  SCREEN. 
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SECTION  PRESSURE  OR AO  COEFFICIENT 


ANGLE  OF  ATTACK,  a  ,  degree* 


Figure  17.  SECTION  PRESSURE  DRAG  COEFFICIENT,^  ,  VERSUS 
ANGLE  OF  ATTACK,  a  ,  IN  UNIFORM  FLOW  WITH 
LOW  SOLIDITY  TURBULENCE  SCREEN. 
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EKPEBIMEMTAL  DATA 
SYMBOL  h( i nchat) 


Figure  18.  NONDIMENSIONAL  SECTION  LIFT,  c.  ,  VERSUS  ANGLE  OF  ATTACK 
NONUNIFORMLY  SHEARED  FLOW. 


Figure  19.  SECTION  LIFT  COEFFICIENT  REFERENCED  TO  MIDCHORD 
FREE-STREAM  DYNAMIC  PRESSURE,  (ct)K  ,  VERSUS  ANGLE 
OF  ATTACK,  or,  NONUNIFORMLY  SHEARED  FLOW. 


SECTION  LIFT  COEFFICIENT  REFERENCED  TO  MIDCHORD  FREE-STREAM  DYNAMIC  PRESSURE,  {€.). 


ANGLE  OF  ATTACK,  CC,  d«grt#» 


Figure  20.  SECTION  LIFT  COEFFICIENT  REFERENCED  TO  HIDCHORD 
FREE-STREAM  DYNAMIC  PRESSURE,  (cj)x  ,  VERSUS  ANGLE 
OF  ATTACK, a,  NONUNIFORMLY  SHEARED  FLOW  (REFERENCE  6). 
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Figure  23.  PRESSURE  COEFFICIENT  REFERENCED  TO  MIDCHORD  FREE-STREAM 
DYNAMIC  PRESSURE,  (CF)n  ,  VERSUS  FRACTION  OF  CHORD,  zfc  , 
FOR  UNIFORM  FLOW  AND  NONUNIFORMLY  SHEARED  FLOW. 
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SECTION  LIFT  COEFFICIENT  REFERENCED  TO  FREE-STREAN  DYNAMIC  PRESSURE  OF  STAGNATION  STREAMLINE,  (cj). 


Figure  24(a).  SECTION  LIFT  COEFFICIENT  REFERENCED  TO  FREE-STREAM 
DYNAMIC  PRESSURE  OF  STAGNATION  STREAMLINE, (cj)s  , 
VERSUS  ANGLE  OF  ATTACK  MEASURED  FROM  ZERO  LIFT,  5  , 
NONUNIFORMLY  SHEARED  FLOW. 
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SECTI 


ANGLE  OF  ATTACK  MEASURED  FROM  ZERO  LIFT,  5r  ,  dagraet 


Figure  2M^(b).  SECTION  LIFT  COEFFICIENT  REFERENCED  TO  MIDCHORD  FREE- 
STREAM  DYNAMIC  PRESSURE,  te,;* ,  VERSUS  ANGLE  OF  ATTACK 
MEASURED  FROM  ZERO  LIFT,  s  ,  NONUNI FORMLY  SHEARED  FLOW. 
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PRESSURE  COEFFICIENT  REFERENCED  TO  FREE-STREAM 
DYNAMIC  PRESSURE  OF  STAGNATION  STREAMLINE,  (Cm). 


AIRFOIL  COORDINATE,  */C  ,  FRACTION  OF  CHORD 


Figure  25.  PRESSURE  COEFFICIENT  REFERENCED  TO  FREE-STREAM 
DYNAMIC  PRESSURE  OF  STAGNATION  STREAMLINE, (CP)S  , 
VERSUS  FRACTION  OF  CHORD,  */c  ,  FOR  a  =  5  DEGREES, 
UNIFORM  FLOW  AND  NONUNIFORMLY  SHEARED  FLOW. 
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PRESSURE  COEFFICIENT  REFERENCED  TO  FREE-STREAM 
DYNAMIC  PRESSURE  OF  STAONATION  STREAMLINE, (Cp)^ 


Figure  26.  PRESSURE  COEFFICIENT  REFERENCED  TO  FREE-STREAM 
DYNAMIC  PRESSURE  OF  STAGNATION  STREAMLINE ,(C,)S  , 
VERSUS  FRACTION  OF  CHORD ,*/c  ,  FOR  &  =  12.5  DEGREES, 
UNIFORM  FLOW  AND  NONUNIFORMLY  SHEARED  FLOW. 
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STREAMLINE  DISPLACEMENT,  AY.  ,  inches 


STREAMLINE  DISPLACEMENT,  AV  ,  inches 


ANGLE  OF  ATTACK,  »  ,  degrees 


Figure  28.  DISPLACEMENT  OF  NONUNIFORMLY  SHEARED  FLOW  AIRFOIL  STAGNATION 
STREAMLINE  RELATIVE  TO  AIRFOIL  MIDCHORD,  h=-2  INCHES. 
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STREAMLINE  DISPLACEMENT,  AV,  ,  inches 


STREAMLINE  DISPLACEMENT,  AYf  ,  inches 


-16-12-8  -4  0  4  8  12  16  20  24 


ANGLE  OF  ATTACK,  QC  ,  degrees 

Figure  30.  DISPLACEMENT  OF  NONUNIFORMLY  SHEARED  FLOW  AIRFOIL  STAGNATION 
STREAMLINE  RELATIVE  TO  AIRFOIL  MIDCHORD,  h  =  +1  INCH. 
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Figure  31.  DISPLACEMENT  OF  NONUNIFORMLY  SHEARED  FLOW  AIRFOIL  STAGNATION 
STREAMLINE  RELATIVE  TO  AIRFOIL  MIDCHORD,  h  =  +3  INCHES. 
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PRESSURE  COEFFICIENT  REFERENCED  TO 
MIDCHORD  FREE-STREAM  DYNAMIC  PRESSURE,^ 
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Figure  32.  COMPARISON  OF  COMPUTED  AND  MEASURED  PRESSURE  COEFFICIENTS, 
(C,)„  VERSUS  FRACTION  OF  CHORD,  x/c  ,  NONUNIFORMLY 
SHEARED  FLOW,  h  =  -2  INCHES,  a  =  7.1  DEGREES. 
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PRESSURE  COEFFICIENT  REFERENCED  TO 
NIDCNORO  FREE-STREAM  DYNAMIC  PRESSURE,  (Cp)^ 


AIRFOIL  COORDINATE,  X/c  ,  FRACTION  OF  CHORD 


Figure  33.  COMPARISON  OF  COMPUTED  AND  MEASURED  PRESSURE  COEFFICIENTS, 
(C, )«,  VERSUS  FRACTION  OF  CHORD,  tic  ,  NONUNIFORMLY 
SHEARED  FLOW,  h  =  -2  INCHES,  a«9.5  DEGREES. 
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PRESSURE  COEFFICIENT  REFERENCED  TO  ,  % 

MIDCHORD  FREE-STREAM  DYNAMIC  PRESSURE,  (C,)K 


Figure  34.  COMPARSION  OF  COMPUTED  AND  MEASURED  PRESSURE  COEFFICIENTS, 
( CP, )„  VERSUS  FRACTION  OF  CHORD,  k/c  ,  NONUNIFORMLY 
SHEARED  FLOW,  h  =  +3  INCHES,  8.2  DEGREES. 
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SECTION  LIFT  COEFFICIENT  REFERENCED  TO  MIDCHORD  FREE-STREAM  DYNAMIC  PRESSURE,  (Cjfo 


-6-4-2  0  2  4  6  8  10  12  14  16  18 


ANGLE  OF  ATTACK,  <Z  ,  degrees 


Figure  35.  COMPARISON  OF  COMPUTED  AND  MEASURED  SECTION  LIFT 
COEFFICIENT  REFERENCED  TO  MIDCHORD  FREE-STREAM 
DYNAMIC  PRESSURE,  NONUNI FORMLY  SHEARED  FLOW. 
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SECTION  QUARTERCHOND  MOMENT  COEFFICIENT 


Figure  36. 


SECTION  LIFT  COEFFICIENT  REFERENCED  TO 
MIDCHORD  FREE-STREAM  DYNAMIC  PRESSURE,  (Cj)K 

COMPARISON  OF  COMPUTED  AND  MEASURED  SECTION  QUARTERCHORD 
MOMENT  COEFFICIENT,  NONUNIFORMLY  SHEARED  FLOW. 
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APPENDIX 


MATHEMATICAL  FORMULATION  AND  IMPLEMENTATION 
OF  SHEARED  FLOW  AERODYNAMICS  COMPUTER  PROGRAM 


FUNDAMENTAL  EQUATIONS 

The  flow  model  on  which  the  analysis  is  based  is  illustrated 
in  Figure  2.  The  wind-tunnel  walls  are  represented  by  vortex  sheets 
of  variable  strengths  (z)  .  where  i  =  1  refers  to  the  lower  wall 
and  4=2  refers  to  the  upper  wall.  Likewise,  the  airfoil  surfaces 
are  represented  by  vortex  sheets  of  variable  strengths  Y^t(x)  f  Yai^» 
and  7alW- 


Point  (a),  Figure  2,  in  the  velocity  profile  marks  the 
boundary,  represented  by  the  function  q(x)  ,  between  different 
constant  values  of  the  vorticity.  In  Figure  2,  the  vorticity  above  this 


boundary  is 


Uj’Vs 

L 


-  u> 


and  below  this  boundary  is 


U. 


00. 


The  functions  ,  y2,  yAL,  fa,  742,and  $  are  the  unknowns  which  are 
to  be  determined. 
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Jj 

The  velocity  components  o,  and  v*  can  be  written 

ouU.y)  *  OL".  +U0  +aA 

*(*.y)  •  *„,  +  *0+  *rA  J 


(A-l) 


where  (  v^.)  are  contributed  by  the  wind-tunnel  wall  singularity 

distribution,  fa^t^arise  from  the  vorticity  u)  in  the  shear  flow, 
and  (oA,vi)  are  the  results  of  the  airfoil  singularity  distribution. 
These  follow  directly  from  application  of  the  Biot-Savart  induction 
law,  and  are  as  follows: 


/ 

2jt 


•  rgtt)d{ 

<Z’t)'  +  (y-L)* 


+ 


(*-()*  + (y+1)2 


(A  -2) 


(z-f)*+(y-L)z 


(x-f)*+(y+L)z 


^ 


-[  ptiV  (y-rf)dff 


m 


<*-<)*+  (y-v)2 


(y-v)*?  1 


(A-3) 


(A-4) 


Symbols  used  only  in  this  appendix  are  defined  as  they  are 
first  introduced. 
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Vw  (X 


aJ  r°°  j  f9(<)  (z-j 

,y)  S  '  2tt  «/„  1  J.L  (*’{)*  + 


9(f)  (z-S)dq 


-f 


(z-€)dy 


<g({)  (x>-4)z  +(<j-qf 


dt 


(A-5) 


**(*.(/)  -  "  — 


i  i  i'*"  -J /+[S('t)Y(‘i-i) 

**  [*-*7)]z+(.ri)‘ 


-MOnl'! 

^-hAUU)]  YuM 

(*-()*  + 

'y-hAU(f)]‘ 

/ 


(A- 6) 


,.rru\.  J-  1  fV*U  ft*  [*-■*(*))]  Yu  dr) 

\J¥m  ■  [w,7)]<  +  (v.7)* - 


fXr  ~Jn-[h'tu(4}V  (*-f>Yydf 
Jz,.  <*-()*+  ['j-h,„(i)Y 

*r 

<*-<)'+  [‘/-hAi(oy 


f. 


(A-7) 


fl/L/y) 

where  h'(x)  5  -^-,  hAu(x)  and  hAL  ( are  y  -  coordinates  of  the 
upper  and  lower  aft  surfaces,  respectively,  of  the  airfoil, 

s  dd(y)  *  stfy)  are  x  -coordinates  of  the  leading  edge 
of  the  airfoil,  ^ L  ,  yM  are  the  y  -coordinates  of  the  lower 
surface  and  upper  surface  points,  respectively,  dividing  the  airfoil 
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leading-edge  vortex  sheet  and  airfoil  upper  and  lower  aft  vortex  sheets 

and  (  Ztg)  are  dummy  variables  of  integration  corresponding  to  (j&. ). 
The  necessity  for  separate  treatment  of  the  airfoil  leading  edge  is  a 
result  of  the  infinite  slope  there  in  terms  of  the  (/-coordinates  of  the 
airfoil  surface. 

BOUNDARY  CONDITIONS 


The  boundary  conditions  which  must  be  satisfied  are 

1.  Velocity  components  normal  to  the  wind-tunnel  walls  are 
zero,  or 

L)  •  0.  (  A-8  ) 

2.  Velocity  components  normal  to  the  airfoil  surface  are  zero  ; 
i.  e.r 

dhA  iS(x,,hA(x))  ,  A  ft% 

‘  (A-9) 

3;  The  boundary  g(z)  is  a  streamline;  i.  e. t 

oLg  _  ■v'(X'tg(x)) 

dx  u(z,gfc))  , 

where  hA  represents  either  hA(J  or  hAL  as  the  case  may  be.  This  condition 
is  satisfied  if 

r9(x)  l 

/  u(x,tj)dg  *  ~  (Ut*  U3)  .  (A- 10) 

-L  * 

Also,  the  Kutta-Joukowski  condition  must  be  satisfied  at  the  airfoil 
trailing  edge;  i.  e.  ,  the  velocity  at  the  trailing  edge  is  finite. 

Equations  (A-8),  (A-9),  and  (A-10)  with  Equation  ( A  - 1 )  result 
in  six  nonlinear  integral  equations  for  the  six  unknown  functions  y, ,  7^, 
7al*  7a,’  7 az  >  and  g  .  It  is  possible  to  obtain  a  solution  only  by 
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recourse  to  a  high-speed,  high-storage-capacity  digital  computer.  The 
nonlinearity  in  the  problem  arises  through  gM  ;  if  g  were  known,  then 
solutions  for  the  y'b  could  be  obtained  by  a  direct  inversion  procedure 
on  the  computer.  The  Kutta- Joukowski  trailing -edge  condition  is 
applied  by  requiring  that  in  the  limit  as  the  airfoil  trailing  edge  is 
approached, the  total  circulation  per  unit  chord  approaches  zero;  i.  e. , 
(see  Figure  2) 


&nv  [yai  (*)  =  0  • 

X  7b  ym 


(A-ll) 


IMPLEMENTATION  OF  DIGITAL  COMPUTER  PROGRAM 


Expressions  for  Velocity  Components 


It  is  assumed  that  the  flow  disturbance  at  the  wind-tunnel  walls 
is  confined  to  0W{ see  Figure  2),  and  that  the  flow  disturbance 

near  the  flow  centerline  (  y  -  0)  is  confined  to  ~DS  <  %  £  0S  .  In  the  far 
stream  (  X  <  -Dw  ,  X  >  0#)  7,  and  y2  are  assumed  constant;  namely, 


y^x)  *  -  0,  , 


for 


fc(x) M  Oz  ,  J 


X/  <  —  Drf ,  x  >  0^ 


For  the  purposes  of  numerical  computation,  the  vortex  sheets 
representing  the  wind-tunnel  walls  and  the  airfoil  surfaces  are  divided 
into  segments,  Lxl  in  length.  The  vortex-sheet  strengths  {  ff  ,  7t  * 
7al  •  7ai  »  and  7a2  in  Figure  2)  are  assumed  to  be  constant  within 
increments  (but,  of  course,  varying  from  segment  to  segment). 
Note  that  because  of  this  assumed  variation  of  vortex  sheet  strength,  v' 
is  singular  at  the  point  where  two  vortex  sheet  segments  of  unequal 
strengths  join. 
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The  region  between  x>--0s  and  k  =  +ds  is  also  divided  into  increments 
For  the  purpose  of  evaluating  the  integrals  of  Equations  (A-2) 
and  (A-3)  numerically,  the  value  of  g(%)  is  assumed  to  be  constant  in 
a  given  increment  Ax±  ,  but  varying  from  increment  to  increment. 


Once  the  vortex  sheets  have  been  reduced  to  segments  of 
constant  singularity  strength  per  unit  length,  the  expressions  for  the 
velocity  components  u  and  xr  (  Equation  (A - 1 ) >  can  be  written 


U(x,i/)  = 


—  cd  LLj 


~  (0L2  +  +  64- ^  ^  &  ) 


(A  - 12) 


and 


viz.ci) 


uj  %r, 


+  v'z  +  v'd  +  f  TTf  *  v'(r  . 


(A- 13) 


Note  that  in  Equations  (A-2)  through  (A-7)  and  in  Equations  (A-12)  and 
(A-l  3), 


m7p(u*  + a*)  •  ^  =  ~hr  (Vz  +  V'3) 

u)  uJ 

Uur  "  ¥rr  U't  ’  ^  "  27F  V'1 

ocA  =  (“**  +  **  +  u<-)  y  v*A  *  (V*  t  i/y  t  w6)  . 
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The  terms  ut  and  ir  are  associated  with  the  various  singularities  as 
follows: 


TERMS 


ASSOCIATED  WITH 


r, 

(A  ^  > 

if* 

» 

. 

*6 

These  terms, 

Distributed  vorticity 
Lower  wind-tunnel  wall  vortex  sheet 
Upper  wind-tunnel  wall  vortex  sheet 
Airfoil  leading-edge  vortex  sheet 
Airfoil  aft  lower -surface  vortex  sheet 
Airfoil  aft  upper  -  surface  vortex  sheet 


individual  vortex  sheet  segments,  become 


Ai* 


ut(x,  y)  *  utD  +  JT  Au-ii 

i  =7 


(A- 14) 


a2(x,,y)  =  -Ut  u2t>  +£  /,iAU'2i 

i*  t 


(A-15) 


u3(x,y)  =  U 2  tciD  *Z  YiiAuii 

i*/ 


(A-16) 


-V, 


(x.y)  -  21  ValL  a 
i’/ 


(A-17) 


us(%,tj)=  yA(i  a  List 


o  *1 


(A-18) 


n2 

*  Z  /A2i 

L*t 


(A- 19) 
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(A-20) 


V,(i t,y)  •  Z 


^.y)  =  -0,  no  +  Z  Yu  &  Vti 

i*  1 

v-sfay)  c  u2  v3p  +  JT  Yzi 

i=t 

* 1 

tfifM  =  22  Yali  &r¥i 
*  ■  / 


*> 

*9<*v)a  22  Yakov's* 
<•/ 

n(*.y)  *Z1  Yah  Av*t 

*•/ 

where 


(z-O)  log 

1  [(z-0)a  +  ylY 

\  [(% -D)*+  (y  +L)*][{z-  D)*+(<j-l)z] 

-(zfD)iog  < 

1 

[(z-tD)*  +  yzY 

,  [(z*D)1  +(y  +  D2]  [(*  +D)Z+  (<J 

-til  i 

+  - 

L )  +  Vyftan  y  'j  -  tan  , 

(¥)] 

+  z(y+L)\ 

ian  (**/)  -*an 

+  l{y-L)  [ 

tan  'iy-l)-*™  J 

(A-21 ) 

(A-22) 

(A-23) 

(A-24) 

(A-25) 


(A-26) 
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u3t>  *  ~7r^-1(^)  ~t(Lr>'i  (-fr r ) 

(*’$**»' 1(z~rfi)  ~(%+D)tour)~1{'^rb)  +  {yty 


(z-D)z  +  y2 
(zi-  D)z+  yz 


(A-27) 

(A-28) 

(A-29) 


Vjo 


1  lo  2  ^  ^ 

2  3  (z+D)z  +  ( y+i)2 


VjfO  3 


_£/  (z-D)z  +(y-L)z 
lU>g  (X  +  D)2  +  (y-L)z 


(A-30) 

(A -31) 


A^/4-  =  (x-z^jAx^log 


(z-zitj&ziy+  (tf<-i)'|i (z-z;  +  j&xif+(y-l)2 


• - 1 

*  >* 

I 

2 

[(*■' 

(*’*i  ’{&*i) 

\l+(vDz\ 

Hy 


-9i)  _ 


.  / *-*.*{&*.  \  , 

Un  v  /  *  <an  l- f* — ; 


,  r  /*  ~  X:  +  j-  A%;  \  ,  (*-*i  -  T  h*L\ 

-Z(y  +  L)  [4an-f[ - - )  •  tan~  ( - yTT - 7_ 


-l(y-L)  \t<zn 


t 


^/z-Zi+jAxj 


y-L 


-  tan 


->ppi) 


A  U  21  ■  ^<a/» 


-/ 


(A-32) 

(A-33) 

(A-34) 


99 


Aa 


1 

i  [v-Ui+A  (*-A)  *  2  ‘  %] 

pd'i)’ 

[*‘4 -.^Y  «/-&)]* 

-A' 


tar» 


-1 


X,  -iL  -li  (  y  -  y-  ) 


~  toil 


-1 


.  »  i+UD* A  ' 

+Ji  (x-^i) - -z - Ay- 

z-A 


(A-35) 


* 


1 

H 

iii  .  i  <  +  (h»uV  .  1 

ta.n~t 

H 

Z’ti+bMiiy-htLi)  +  2 

i '*(&!)* 

trhALi~hALi 

—  tan 


-t 


*-*:  +  hui(4’hALi)-  f  +  (*AL^  a»; 

y-hALl  -hALi(*-*i) 


\x-X;+hAa(y- hAU)  +  2  d*;  ] 

M 

y~^Au  ~^au 

Z'^i^^ALi  (y-htii)-  £'  A*,I 

2 

•f 

[y-^Li-hAufx-x^ 

(A-36) 


A(X^-  is  the  same  expression  as  A^.with  hAl-  ,  hALl_  replaced  by 
h AUt  ">  ^ AUt  • 


A  Vfi  =  (x  -  *i  +  j  A  X  ■  j  tun 


+  l(y-Si)io9 


+  (*-*;  *-j**i)2 

(y-9i)z  *  (*-*>- U*i)z 


(A-37) 
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(z-Zi  +  fAxi)  +  (y+L)' 

A\ru  =  ylog - 7 - ^ - 

(*-*£  -jAZi)  +  (y+0 


(A-38) 


/  (x-Zi+itoif+te-1)' 


(A-39) 


lt  +  U/)z 


~  tcLn 


,  y-i/c+4'(z-A)  +■  1+(£—' 

x-Jt  -j/fy-K) 


y- &  +A'(*  -A)  -  1%2- 

*-4 


2 

+ 

I  1 

N 

i 

1  N 

[</-</i+4(z-4)-  2 

2 

+■ 

[x-4-j/(  «-*>]* 

(A-40) 


^  2/^  —  ^  ^  .  ■ 
SC  2  [ *  *.  *  /»4U  ( A X  J  +  </-  *Ut*  * *  * ’*>] 


/  .  ^  (y  ~  ^ ALL  )  *  2~~~ 

~  - — — - 

y-^a  -  (*'*J 


(cf-hAU)-  UShuL  axl 

i  ■  - 

y-^ALL  -  h'ALi  (*-*;) 


AiriL  is  the  same  expression  as  with  ,  /j4it;  replaced  by 


(A-41) 


^Aui  *  ^Aui  • 
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In  these  expressions  (Equations  (A-14)  through  (A-41)),  a/,  is  the 
number  of  increments  into  which  the  region  <»  6DS  is  subdivided, 
a*2  is  the  number  of  segments  into  which  each  of  the  wind-tunnel  wall 
vortex  sheets  between  -Dw  *  »  <  D#  is  subdivided,  ^  is  the  singularity 
strength  of  the  i**vortex  sheet  segment,  Ni  is  the  number  of  segments 
into  which  the  airfoil  leading-edge  vortex  sheet  is  subdivided,  is  the 
number  of  segments  into  which  the  vortex  sheets  representing  the  aft 
upper  and  lower  surfaces  of  the  airfoil  are  each  subdivided,  (  ,  g-  ) 

are  the  coordinates  of  the  midpoint  of  the  segment  of  yfc)  between 
-  0S  <  x>  <  Ps  •  ( 25t ,  t  L  )  are  coordinates  of  the  midpoint  of  the  i*h 
segment  of  wind-tunnel  wall  vortex  sheet,  (  z;  ,  )  and  (  zt  ,  hAn  ) 

are  the  coordinates  of  the  midpoint  of  the  segments  of  airfoil  upper 
and  lower  vortex  sheets,  respectively,  ( Zi ,  y;  )  are  the  coordinates  of 
the  midpoint  of  the  ith  segment  of  airfoil  leading-edge  vortex  sheet,  and 
{A*;,  Jy;  )  are  the  projections  on  the  x  and  </  axes  of  the  segment  length. 
The  prime  denotes  differentiation  ;  e.  g.  , 


/  _  dh  _ 
at* 


The  principal  values  of  the  inverse  tangents  is  to  be  understood, 
and  logarithms  are  natural  logarithms.  The  general  computational 


procedure  adopted  is  as  follows: 


1,  An  initial  shape  of  the  dividing  streamline,  ql0)(x,),is  assumed. 


2.  The  boundary  conditions  that  the  flow  is  tangential  at  the 
wind-tunnel  walls  and  the  airfoil  surface,  together  with 
the  approximation  to  g(v) ,  are  used  to  determine  the 
singularity  strength  distribution  of  the  vortex  sheets 
representing  the  walls  and  airfoil  surfaces.  The  boundary 
conditions  at  the  wind-tunnel  walls  and  at  the  airfoil  sur¬ 
face  are  satisfied  at  discrete  points,  x  •  >  corresponding 
to  the  segments  Ay^  of  vortex  sheets  along  these  surfaces. 
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3.  The  second  approximation,  is  determined  on  the 

basis  of  the  mass  flow  continuity  requirement  based  on 
Equation  (A  -10), 

4.  The  iteration  proceeds  for  the  adjusted  boundary  by 
going  back  to  Step  2  with  the  new  approximation  to  g(z). 

Boundary  Conditions  -  "Smoothing"  Technique 

When  the  velocity  components  u  and  xr,  Equations  (A- 12)  through 
(A-41),  are  substituted  into  the  equations  for  the  boundary  conditions, 
Equations  (A-8)  and  (A-9),  and  including  Equation  (A-ll),  there 
results  *  Ay  *  2rtz  +  /  linear  equations  in  terms  of  Zr^2  *  H 
unknown  Yi  1  s.  This,  of  course,  is  an  overdeterrnined  set  of  equations. 
(As  noted  in  the  body  of  this  report,  this  overdeterminacy  is  very 
likely  a  result  of  the  particular  assumptions  inherent  in  the  method  of 
numerical  analysis  adopted;  in  particular,  constant  singularity  strength 
along  small  segments  of  the  vortex  sheets.  )  The  "solution”  of  this  set 
of  equations  is  effected  by  means  of  a  smoothing  process,  based  on  the 
results  presented  in  References  10  and  11.  This  approach  was  resorted 
to  because  of  the  unsatisfactory  behavior  of  solutions  obtained  by  direct 
inversion  of  the  system  of  equations  (see  Figures  3  and  4).  The  set  of 
equations  can  be  written 

H  t4  f) 

L  Atj  Yn  +  II  Bcj  Yzl  +  ZT  c-  y3i  =  Dj  .  (A-42) 

i  ~  /  i,-t  i  *  f 

(The  Yji  are  values  of  y  on  the  airfoil,  and  H  -  A/f  *  2a/2  . )  It  is 
assumed  that  the  ^  are  in  error  by  some  small  amount,  £j  .  At  the 
same  time,  it  is  known  from  experiment  that  the  //;  ,  ,  and  fsi  a*e 

smoothly  varying  qxiantities.  To  achieve  this  in  the  mathematical 
solution,  it  is  proposed  to  minimize  second  differences  of  the  Yi  .  Thus 
it  is  desired  that 
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Ml  -  /  Ml  -  T 

ZZ!  (ifid-i)  ~  2  Yu*  Yui+  /))  *  ZZ  (  Yui-()~  2  yZi  +  /zu+f)) 
c*  2 


Mi-1 

Z 

i  «•  2 


n 

+  y' (iidti)  *-f3i  +  /j(n-i))z 
i  */ 


be  a  minimum,  subject  to 

(2* li+n+1) 

E  *' 


-  constant. 


The  Yji  are  continuous  around  the  airfoil  surface; 
is  adopted  that  y30  s  y3n ,  ^  ^ to  insure  this. 


hence,  the  convention 


The  development  follows  closely  that  outlined  in  Reference  11, 
and  details  are  omitted  here.  The  results  can  be  stated  simply  in 
matrix  form.  If  Equation  (A-42)  is  written  in  matrix  form 

Oily]  *  L°J.  <A-43> 

then  the  desired  solution  is  obtained  from 

[A*A  +KH][r]  =  M[P]  (A-44) 


where  the  matrix  [  A*  ]  is  the  transpose  of  the  matrix  [A  ],  K  is  an 
arbitrary  constant,  and  the  matrix  [//]  is 


0 

\  0 

M- 

0 

*i 

!  (9 

1 

0 

L. 

0 

|  Hz  _ 

a  ( ltft  *  n  )  by  (  2Mz+n)  matrix,  where 


(A-45) 
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an  by  matrix,  and 


an  n  by  n  matrix. 
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If  K  -  0,  the  solution  obtained  represents  a  "least  squares" 
optimum  solution  with  no  smoothing.  As  K  is  increased  from  a  zero 
value,  more  and  more  smoothing  is  evident  in  the  solutions  obtained, 
but  in  general  the  sum  becomes  larger  also.  The  uniform  flow 
solution  shown  in  Figure  4  was  obtained  with  K  =  0.  1 ;  the  sheared  flow 
solutions  (Figures  32  through  34)  were  computed  with  K  =  0.  3.  An 
indication  of  the  order  of  magnitude  of  the  cL  is  afforded  for  any 
particular  solution  by  noting  the  airfoil  trailing  edge  values  of  7^  and  yA, 
(Equation  (A-ll));  in  a  typical  calculation,  they  differed  in  magnitude 
by  approximately  2  percent  and  £~  0.003. 


Iterative  Technique 

Given  the  input  values  of  f)  for  the  iteration  and  once  the 
system  of  equations  for  the  ^  has  been  solved  based  on  tno 

velocity  components  u  and  v-  can  be  computed  from  Equations  (.  1-J2) 
through  (A -41)  for  any  value  of  x  and  y  in  the  flow.  The  techniques  used 
in  the  computer  program  to  obtain  input  values  g/J)ior  the  (j *■*)**  iteration 
was  as  follows.  At  the  je4-  centered  on  segments  corresponding  to  the 
constant  ,  the  integral  equation  (mass  flow  continuity) 


+(j-V 

=  r(U,*U3)  (A -48) 

(*.tM  2 

is  solved  numerically  for  The  term  (ir-ug')  =  0  if  the  flow 

normal  to  the  curve  ^*f.)is  zero.  The  input  for  the  ( J+t  )  iteration 
is  then 


J  ucLy  +  (V”U.g-) 
-L 


(A-49) 
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Initially,  calculations  were  made  using 

d*  -  S'0"’ 

but  these  calculations  failed  to  converge. 


Once  the  calculation  has  resulted  in  a  converged  solution,  the 
stagnation  streamline  dynamic  pressure  is  determined  by  computing 
the  mass  flow  between  the  airfoil  surface  and  the  wind-tunnel  wall;  i.  e.  , 

phL(*>) 

J  u(x>,  y)d<f  . 

-L 

Finally,  the  airfoil  stagnation  streamline  dynamic  pressure  is  p[u(y^j2 
where 


r4*  rh^<xi 

j  U(y)cLy  =  /  u(x,tj)d<4 

- l  -  L 


♦ 
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